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Preface 


View a promotional video for this course on YouTube 


These are my lecture notes for my online Coursera course, Differential Equations for Engineers. I 
have divided these notes into chapters called Lectures, with each Lecture corresponding to a video on 
Coursera. I have also uploaded all my Coursera videos to YouTube, and links are placed at the top of 
each Lecture. 

There are problems at the end of each lecture chapter and I have tried to choose problems that 
exemplify the main idea of the lecture. Students taking a formal university course in differential 
equations will usually be assigned many more additional problems, but here I follow the philosophy 
that less is more. I give enough problems for students to solidify their understanding of the material, 
but not too many problems that students feel overwhelmed and drop out. 1 do encourage students to 
attempt the given problems, but if they get stuck, full solutions can be found in the Appendix. 

There are also additional problems at the end of coherent sections that are given as practice quizzes 
on the Coursera platform. Again, student should attempt these quizzes on the platform, but if a 
student has trouble obtaining a correct answer, full solutions are also found in the Appendix. 

Students who take this course are expected to know single-variable differential and integral calcu- 
lus. Some knowledge of complex numbers, matrix algebra and vector calculus is required for parts of 


this course. Students missing this latter knowledge can find the necessary material in the Appendix. 


JEFFREY R. CHASNOV 
Hong Kong 
January 2019 
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Lecture 1 


Introduction to differential equations 


View this lecture on YouTube 


A differential equation is an equation for a function containing derivatives of that function. For exam- 
ple, the differential equations for an RLC circuit, a pendulum, and a diffusing dye are given by 


dq pa, 1 e 
La H Rz tod Eycoswt, (RLC circuit equation) 
Po do . 
mi qa? cl gp T "89 0 = Focoswt, (pendulum equation) 
du du u Pu en 
ay D ( e =) . (diffusion equation) 


These are second-order differential equations, categorized according to the highest order derivative. 

The RLC circuit equation (and pendulum equation) is an ordinary differential equation, or ode, and 
the diffusion equation is a partial differential equation, or pde. An ode is an equation for a function of 
a single variable and a pde for a function of more than one variable. A pde is theoretically equivalent 
to an infinite number of odes, and numerical solution of nonlinear pdes may require supercomputer 
resources. 

The RLC circuit and the diffusion equation are linear and the pendulum equation is nonlinear. In 
a linear differential equation, the unknown function and its derivatives appear as a linear polynomial. 
For instance, the general linear third-order ode, where y = y(x) and primes denote derivatives with 
respect to x, is given by 


m 


ay” + a(x)y" + ay(x)y! + ao(x)y = b(x), 
where the a and b coefficients can be any function of x. The pendulum equation is nonlinear because 
of the term sin@, where O = 0(t) is the unknown function. Making the small angle approximation, 
sin 6 ~ 6, the pendulum equation becomes linear. 
The simplest type of ode can be solved by integration. For example, a mass such as Newton’s 
apocryphal apple, falls downward with constant acceleration, and satisfies the differential equation 


d?x 


de 


With initial conditions specifying the initial height of the mass xy and its initial velocity ug, the solution 
obtained by straightforward integration is given by the well-known high school physics equation 


1 
x(t) = xo + uot — aah 


1 


LECTURE 1. INTRODUCTION TO DIFFERENTIAL EQUATIONS 


Problems for Lecture 1 


1. Classify the following differential equations by order, ordinary or partial, and linear or nonlinear. 


dey dy 
kae ae 
1d 2dp\ _ $ 
ER 6 A 
ö du l yu _ ou 
ot ox dx? 
d?x dx 
d) aa FO rer 0 


u (5 u 2) 
e) | | 


Solutions to the Problems 


Week I 


First-order differential equations 


In this week’s lectures, we discuss first-order differential equations. We begin by explaining the 
Euler method, which is a simple numerical method for solving an ode. Not all first-order differential 
equations have an analytical solution, so it is useful to understand a basic numerical method. Then the 
analytical solution methods for separable and linear equations are explained. We follow the theory and 
some simple examples with three real-world applications of first-order equations and their solution: 


compound interest, terminal velocity of a falling mass, and the resistor-capacitor electrical circuit. 


Lecture 2 


Euler method 


View this lecture on YouTube 


Y, = Yy+AX (Xp Y q) 


| \ 


slope = X Yo) 


0 1 


The ode dy/dx = f(x,y) with initial condition y(xo) = yo is integrated to x = x, using the Euler method. 


Sometimes there is no analytical solution to a first-order differential equation and a numerical solution 
must be sought. The first-order differential equation dy/dx = f (x,y) with initial condition y(xo) = Yo 
provides the slope f(xo,yo) of the tangent line to the solution curve y = y(x) at the point (xo, yo). 
With a small step size Ax = x] — xp, the initial condition (xo, yo) can be marched forward to (x1,y1) 
along the tangent line using Euler’s method (see figure): 


yı = Yo + Axf (xo, y0). 


This solution (x1,y1) then becomes the new initial condition and is marched forward to (x», y2) along 
a newly determined tangent line with slope given by f(x1,y1). For small enough Ax, the numerical 
solution converges to the unique solution, when such a solution exists. 

There are better numerical methods than the Euler method, but the basic principle of marching the 


solution forward remains the same. 


8 LECTURE 2. EULER METHOD 


Problems for Lecture 2 
1. The Euler method for solving the differential equation dy/dx = f(x,y) can be written in the form 
ky =Axf(Xn,Yn), Yny = Yn + kı, 
and is called a first-order Runge-Kutta method. Second-order Runge-Kutta methods have the form 
ky = Axf(Xn,Yn), k2 = Axf (Xn +0Ax,Yn + ki), Yn+1 = Yn + akı + bkz, 


where a +b = 1 and ab = Bb = 1/2. Write down the second-order Runge-Kutta methods with a = b 
and with a = 0. These methods are called the modified Euler method (or predictor-corrector method) 
and the midpoint method, respectively. 


Solutions to the Problems 


Lecture 3 


Separable first-order equations 


View this lecture on YouTube 


A first-order ode is separable if it can be written in the form 


gW) = f(x), yo) = yo, 


where the function g(y) is independent of x and f(x) is independent of y. Integration from xo to x 


results in 


[ suey @ax = [far 


xo 
The integral on the left can be transformed by substituting u = y(x), du = y'(x)dx, and changing the 
lower and upper limits of integration to y(x9) = yo and y(x) = y. Therefore, 


[cod = [fax 


which can often yield an analytical expression for y = y(x) if the integrals can be done and the 
resulting algebraic equation can be solved for y. 
A simpler procedure that yields the same result is to treat dy/dx as a fraction. Multiplying the 


differential equation by dx results directly in 


g(y)dy = f(x)dx, 


which is what we call a separated equation with a function of y times dy on one side, and a function 
of x times dx on the other side. This separated equation can then be integrated directly over y and x. 


10 LECTURE 3. SEPARABLE FIRST-ORDER EQUATIONS 


Problems for Lecture 3 


1. Put the following equation in separated form. Do not integrate. 


a Y - x?y — 4y 
dx x+4 


b) E = sec(y)e* ¥(1 + x) 


dy _ xy 
o) dx (x+1)(y+1) 


do 
d) — +sin@ =0 


Solutions to the Problems 


Lecture 4 


Separable first-order equation: example 


View this lecture on YouTube 
Example: Solve y +y?sinx=0, y(0) =1. 


We first manipulate the differential equation to the form 


and then treat dy/dx as if it was a fraction to separate variables: 


d 
Mn —sin x dx. 


y? 


We then integrate the right side from x equals 0 to x and the left side from y equals 1 to y. We obtain 


y dy Pie 
I ==- sin x dx. 
1 y 0 


Integrating, we have 


y x 
—-| =cosx| , 
1 0 
or 
1— -= = cosx- 1. 
Solving for y, we obtain the solution 
1 
Y= 3 cosx 


11 


12 LECTURE 4. SEPARABLE FIRST-ORDER EQUATION: EXAMPLE 


Problems for Lecture 4 

1. Solve the following separable first-order equations. 
a) dy/dx = 4x,/y, with y(0) = 1. 
b) dx/dt = x(1 — x), with x(0) = xy and 0 < xp < 1. 


Solutions to the Problems 


Practice quiz: Separable first-order 


odes 


1. Find the solution of y” = ,/Xy, with initial value y(1) = 0. 
2. Find the solution of y? — xy’ = 0, with initial value y(1) = 1. 
3. Find the solution of y’ + (sin x)y = 0, with initial value y(rr/2) = 1. 


Solutions to the Practice quiz 
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14 


LECTURE 4. SEPARABLE FIRST-ORDER EQUATION: EXAMPLE 


Lecture 5 


Linear first-order equations 


View this lecture on YouTube 


A linear first-order differential equation with initial condition can be written in standard form as 


Ys p(xy= s(x), yl) = yo: 6.1) 


All linear first-order differential equations can be integrated using an integrating factor u. We multiply 


the differential equation by the yet unknown function u = u(x) to obtain 


nx) | + play] = usa); 


and then try to determine u(x) so that 


The unknown function u(x) is called an integrating factor because the resulting differential equation, 


Kamen) = u(x)g(x), can be directly integrated. Using y(xo) = yo and choosing u(xo) = 1, we have 


u(x)y — yo = [ u(x)g(x)dx; or after solving for y = y(x), 
1 x 
y) = (vo + f nagoa): 62) 


d 
To determine u(x), we differentiate and expand Em [n(x)y] = u(x)g(x) to yield 


dy _ dy dy 
FR + puy = q 


and upon canceling terms, we obtain the differential equation 


du 
= p(x), na) =1 


This equation is separable and can be easily integrated to obtain 


u(x) = exp (f pax), (5.3) 


0 


Equations (5.2) and (5.3) together solve the first-order linear equation given by (5.1). 


15 


16 LECTURE 5. LINEAR FIRST-ORDER EQUATIONS 


Problems for Lecture 5 


1. Write the following linear equations in standard form. 


2. Consider the nonlinear differential equation dx/dt = x(1 — x). By defining z = 1/x, show that the 
resulting differential equation for z is linear. 


Solutions to the Problems 


Lecture 6 


Linear first-order equation: example 


View this lecture on YouTube 
dy ee 
Example: Solve de +2y =e *, with y(0) = 3/4. 


Note that this equation is not separable. With p(x) = 2 and g(x) = e *, we have 


u(x) = exp UN 2d) 


and 


which can be simplified to 


17 


18 


Problems for Lecture 6 


1. Solve the following linear odes. 


a) ru y(0)=-1 
») Y =2x(1-y), y(0)=0. 


Solutions to the Problems 


LECTURE 6. LINEAR FIRST-ORDER EQUATION: EXAMPLE 


Practice quiz: Linear first-order odes 
1. Find the solution of (1 + x?)y’ + 2xy = 2x, with initial value y(0) = 0. 

2. Find the solution of x?y’ = 1 — 2xy, with initial value y(1) = 2. 

3. Find the solution of y + Ay = a, A > 0, with initial value y(0) = 0. 


Solutions to the Practice quiz 


19 


20 


LECTURE 6. LINEAR FIRST-ORDER EQUATION: EXAMPLE 


Lecture 7 


Application: compound interest 


View this lecture on YouTube 


The compound interest equation arises in many different engineering problems, and we consider 
it here as it relates to an investment. Let S(t) be the value of the investment at time t, and let r be the 
annual interest rate compounded after every time interval At. Let k be the annual deposit amount (a 
negative value indicates a withdrawal), and suppose that a fixed amount is deposited after every time 
interval At. The value of the investment at the time t + At is then given by 


S(t + At) = S(t) + (rAt)S(t) + kAt, 7.1) 


where at the end of the time interval At, rAtS(t) is the amount of interest credited and kAt is the 
amount of money deposited. 
Rearranging the terms of (7.1) to exhibit what will soon become a derivative, we have 


S(t + At) — S(t) 
At 


=rS(t)+k. 
The equation for continuous compounding of interest and continuous deposits is obtained by taking 


the limit At > 0. The resulting differential equation is 


ds 
stk (7.2) 


which can solved with the initial condition S(0) = Sọ, where Sy is the initial capital. The differential 
equation is linear and the standard form is dS/dt — rS = k, so that the integrating factor is given by 


eS", 


The solution is therefore 


t 
S(t) = e" (so +/ keria) 
0 


— Soe” ett (1 a 7 


where the first term on the right-hand side comes from the initial invested capital, and the second 
term comes from the deposits (or withdrawals). Evidently, compounding results in the exponential 


growth of an investment. 
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22 LECTURE 7. APPLICATION: COMPOUND INTEREST 


Problems for Lecture 7 


1. Suppose a 25 year-old plans to set aside a fixed amount each year until retirement at age 65. How 
much must he/she save each year to have $1,000,000 at retirement? Assume a 6% annual return on the 
saved money compounded continuously. Out of the $1,000,000, approximately how much was saved, 
and how much was earned on the investment? 


Solutions to the Problems 


Lecture 8 


Application: terminal velocity 


View this lecture on YouTube 


Using Newton’s law, we model a mass m free falling under gravity but with air resistance. We as- 
sume that the force of air resistance is proportional to the speed of the mass and opposes the direction 
of motion. We define the x-axis to point in the upward direction, opposite the force of gravity. Near 
the surface of the Earth, the force of gravity is approximately constant and is given by —mg, with 
g = 9.8m/s? the usual gravitational acceleration. The force of air resistance is modeled by —kv, where 
v is the vertical velocity of the mass and k is a positive constant. When the mass is falling, v < 0 and 
the force of air resistance is positive, pointing upward and opposing the motion. The total force on 
the mass is therefore given by F = —mg — kv. With F = ma and a = dv/dt, we obtain the differential 
equation 


dv 
ma, = —mg — kv. (8.1) 


The terminal velocity vo of the mass is defined as the asymptotic velocity after air resistance balances 
the gravitational force. When the mass is at terminal velocity, dv/dt = 0 so that 


ae 
vo = -7 (8.2) 


The approach to the terminal velocity of a mass initially at rest is obtained by solving (8.1) with initial 
condition v(0) = 0. The equation is linear and the standard form is dv/dt + (k/m)v = —g, so that the 


integrating factor is 
p(t) = elm, 


and the solution to the differential equation is 


t 
v(t) Sia | Im dt 
0 
__Mg (1 _ „kt/m 
= i (1 e ). 


Therefore, 0 = Vo (1 — et! e and v approaches væ as the exponential term decays to zero. 
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24 LECTURE 8. APPLICATION: TERMINAL VELOCITY 


Problems for Lecture 8 


1. A male skydiver of mass m = 100kg with his parachute closed may attain a terminal velocity of 
200 km/hr. How long does it take him to attain one-half his terminal velocity, and how long to attain 
95%? 


Solutions to the Problems 


Lecture 9 


Application: RC circuit 


View this lecture on YouTube 


RC circuit diagram. 


Consider a resister R and a capacitor C connected in series as shown in the figure. A battery connects 
to this circuit by a switch, stepping up the voltage by €. Initially, there is no charge on the capacitor. 
When the switch is thrown to a, the battery connects and the capacitor charges. When the switch 
is thrown to b, the battery disconnects and the capacitor discharges, with energy dissipated in the 


resister. Here, we determine the voltage drop across the capacitor during charging and discharging. 


The equations for the voltage drops across a capacitor and a resister are given by 
Vo =q/C, Vr=ik, (9.1) 
where C is the capacitance and R is the resistance. The charge q and the current i are related by 


. dg 
= 9.2 
er (9.2) 
Kirchhoff’s voltage law states that the emf E in any closed loop is equal to the sum of the voltage 
drops in that loop. Applying Kirchhoff’s voltage law when the switch is thrown to a results in 


VR+Ve=E. (9.3) 


Using (9.1) and (9.2), the voltage drop across the resister can be written in terms of the voltage drop 


across the capacitor as 
dVc 
Vr = RC, 
ds dt 


25 


26 LECTURE 9. APPLICATION: RC CIRCUIT 


and (9.3) can be rewritten to yield the linear first-order differential equation for Vc given by 


dV, 
a + Vc/RC = €/RC, (9.4) 


with initial condition V¿(0) = 0. 
The integrating factor for this equation is 
u(t) = e RT, 
and (9.4) integrates to 
t 
Ve(t) = e7t/RO f (E/RC)e/ PCat, 
0 


with solution 
Ve(t) =E (1 = eve 


The voltage starts at zero and rises exponentially to €, with characteristic time scale given by RC. 
When the switch is thrown to b, application of Kirchhoff’s voltage law results in 


Vr + Ve =0, 


with corresponding differential equation 


dVc 
Ge + Ve/RC = 0. 


Here, we assume that the capacitance is initially fully charged so that Vc(0) = €. The solution, then, 
during the discharge phase is given by 


Vc(t) = Ee tf RC, 


The voltage starts at € and decays exponentially to zero, again with characteristic time scale given by 
RC. 


Problems for Lecture 9 
1. Determine the current in an RC circuit during charging and discharging. 


Solutions to the Problems 


27 
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LECTURE 9. APPLICATION: RC CIRCUIT 


Practice quiz: Applications 


1. Suppose a 25 year-old plans to set aside a fixed amount each year until retirement at age 65. 
Approximately, how much total must he/she save to have $1,000,000 at retirement? Assume a 10% 


annual return on the saved money compounded continuously. 


2. A male skydiver of mass 100 kg with his parachute closed may attain a terminal velocity of 200 
km/hr. How long does it take him to attain a velocity of 150 km/hr? 


3. An RC circuit consists of a resistor (3000 ©) and a capacitor (0.001 F), where Q is the ohm with 
units kg - m2.s573.A72, Fis the Faraday with units s*.A?2.m”2. kg, and where kg is kilogram, m 
is meters, s is seconds, and A is amps. How long does it take for the capacitor to charge to 95% of the 


battery voltage? 


Solutions to the Practice quiz 
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LECTURE 9. APPLICATION: RC CIRCUIT 


Week II 


Second-order differential equations 


33 


This week’s lectures are on second-order differential equations. We begin by generalizing the 
Euler numerical method to a second-order equation, to show how solutions can be obtained even for 
equations that have no analytical solution. We then develop two theoretical concepts used for linear 
equations: the principle of superposition, and the Wronskian. Armed with these concepts, we can then 
find analytical solutions to a homogeneous second-order ode with constant coefficients. We make use 
of an exponential ansatz, and convert the differential equation to a second-order polynomial equation 
called the characteristic equation of the ode. The characteristic equation may have real or complex 
roots and we discuss the solutions for these different cases. We then consider the inhomogeneous 
differential equation, and the phenomena of resonance, where the forcing frequency is equal to the 
natural frequency of the oscillator. Finally, some interesting and important applications are discussed. 


34 


Lecture 10 


Euler method for higher-order odes 


View this lecture on YouTube 


In practice, most higher-order odes are solved numerically. For example, consider the general second- 
order ode given by 
ž = f (t,x, 2). 


Here, we adopt the widely used physics notation x = dx/dt and & = d?x/dt?. The dot notation is used 
to represent only time derivatives and can be applied to any dependent variable that is a function of 
time. 

To solve numerically, we convert the second-order ode to a pair of first-order odes. Define u = x, 


and write the first-order system as 


x=u, (10.1) 
ú = f(t,x,u). (10.2) 


The first ode, (10.1), gives the slope of the tangent line to the curve x = x(t); the second ode, (10.2), 
gives the slope of the tangent line to the curve u = u(t). Beginning at the initial values (x,u) = (xo, uo) 
at the time t = to, we move along the tangent lines to determine x; = x(to + At) and u1 = u (tọ + At): 


xy = Xo + Atuo, 
uy = uo + Atf (to, Xo, ug). 


The values x1 and u4 at the time ft; = tọ + At are then used as new initial values to march the solution 
forward to time t2 = tı + At. When a unique solution of the ode exists, the numerical solution 


converges to this unique solution as At — 0. 
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36 LECTURE 10. EULER METHOD FOR HIGHER-ORDER ODES 


Problems for Lecture 10 


1. Write the second-order ode, 6 + 6/q + sin0 = f cos wt, as a system of two first-order odes. 


2. Write down the second-order Runge-Kutta modified Euler method (predictor-corrector method) for 
the following system of two first-order odes: 


i=flt,x,y) y=8(t,x,y). 


Solutions to the Problems 


Lecture 11 


The principle of superposition 


View this lecture on YouTube 


Consider the homogeneous linear second-order ode given by 
¥+ p(t)x +q(t)x = 0; 


and suppose that x = Xı(t) and x = X(t) are solutions. We consider a linear combination of Xı and 
X» by letting 

x = (1X1 (t) + c2X2(t), 
with cı and c3 constants. The superpositon principle states that x is also a solution to the homogeneous 


ode. To prove this, we compute 


E+ pt + qx = cıXı + coXo + p (cıXı + c2X2) + q (c1X1 + c2X2) 
cy (Xy +pX1+qX1) +02 (Xo + pX2 +qX>) 
=(1x0+c2Xx0 

=0, 


since X] and Xz were assumed to be solutions of the homogeneous ode. We have therefore shown that 


any linear combination of solutions to the homogeneous linear second-order ode is also a solution. 
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38 LECTURE 11. THE PRINCIPLE OF SUPERPOSITION 


Problems for Lecture 11 


1. Consider the inhomogeneous linear second-order ode given by 


¥+ plt)x + q(t)x = g(t); 


and suppose that x = x„(t) is a solution of the homogeneous equation and x = xp(t) is a solution 
of the inhomogeneous equation. Prove that x = x„(t) + xp(t) is a solution of the inhomogeneous 
equation. 


Solutions to the Problems 


Lecture 12 


The Wronskian 


View this lecture on YouTube 


Suppose that we have determined that x = Xı(t) and x = X>(t) are solutions to 

X+ p(t)x +q(t)x =0, 
and that we now attempt to write the general solution to the ode as 

x =C1X1(t) + c2X2(t). 
We must then ask whether this general solution can satisfy the two initial conditions given by 

x(to) = xo, X(to) = uo. 
Applying these initial conditions to our proposed general solution, we obtain 
c1X1(to) + c2X2(to) = xo, c1X1(to) + ¢2X2(to) = uo, 

which is a system of two linear equations for cı and c2. A unique solution exists provided 


X1 (to) os = X1(to)X2(to) — X1 (to) Xa(to) #0, 


where the determinant given by W is called the Wronskian. 


For example, with w # 0, the two solutions Xı(t) = sinwt and X2(t) = cos wt have a nonzero 


Wronskian for all t since 


W = (sin wt) (-wsin wt) — (w cos wt) (cos wt) = —w. 


When the Wronskian is not equal to zero, we say that the two solutions X;(t) and X(t) are linearly 


independent. The concept of linear independence is borrowed from linear algebra, and indeed, the 


set of all functions that satisfy a second-order linear homogeneous ode can be shown to form a two- 


dimensional vector space. 
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40 LECTURE 12. THE WRONSKIAN 


Problems for Lecture 12 


1. Show that X,(t) = exp (at) and X>(t) = exp (Pt) have a nonzero Wronskian for all t provided 
ar PB. 


Solutions to the Problems 


Lecture 13 


Homogeneous second-order ode with 


constant coefficients 


View this lecture on YouTube 


We now study solutions of the homogeneous constant-coefficient ode, written as 
ax +bx+cx=0, 


with a, b, and c constants. Our solution method finds two linearly independent solutions, multiplies 
each of these solutions by a constant, and adds them. The two free constants are then determined 
from initial values on x and x. 

Because of the differential properties of the exponential function, a natural ansatz, or educated 
guess, for the form of the solution is x = e't, where r is a constant to be determined. Successive 
2ort 


differentiation results in x = re’! and & = 1?e"*, and substitution into the ode yields 


are! + bre’ + cet = 0. 


Our choice of exponential function is now rewarded by the explicit cancelation of e”?. The result is a 
quadratic equation for the unknown constant r: 


ar +br+c=0. 


Our ansatz has thus converted a differential equation for x into an algebraic equation for r. This 
algebraic equation is called the characteristic equation of the ode. Using the quadratic formula, the two 


solutions of the characteristic equation are given by 


a —b + Vv b2 — 4ac 
E 2a " 


There are three cases to consider: (1) if b? — 4ac > 0, the roots are distinct and real ; (2) if b? — 4ac < 0, 
the roots are distinct and complex-conjugates; and (3) if b? — 4ac = 0, the roots are repeated. We will 


consider these three cases in turn. 
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42 LECTURE 13. HOMOGENEOUS SECOND-ORDER ODE WITH CONSTANT COEFFICIENTS 


Problems for Lecture 13 

1. For the following differential equations, determine the roots of the characteristic equation. 
a) x-x=0; 
b) i+x=0; 
c) x-2x+x=0. 


Solutions to the Problems 


Lecture 14 


Case 1: distinct real roots 


View this lecture on YouTube 


When the roots of the characteristic equation are distinct and real, then the general solution to the 


second-order ode can be written as a linear superposition of the two solutions e” 1% and e”2f; that is, 
x(t) = ce"! + eze”. 


The unknown constants cı and cz can then be determined by the given initial conditions x(to) = xo 


and x(fo) = uo. 

Example: Solve ¥ + 5x + 6x = 0 with x(0) = 2, x(0) = 3. 

We take as our ansatz x = e and obtain the characteristic equation 
r? +5r+6=0, 


which factors to 
(r+3)(r +2) =0. 


The general solution to the ode is thus 
x(t) =ac + ce”. 
The solution for x obtained by differentiation is 
x(t) = ut, 
Use of the initial conditions then results in two equations for the two unknown constant cı and cp, 
a+n=2, —2cy — 3c2 = 3, 


with solution cı = 9 and ca = —7. Therefore, the unique solution that satisfies both the ode and the 


initial conditions is 


x(t) = ze 


= 9e 2! ( — Ze) : 
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44 LECTURE 14. CASE 1: DISTINCT REAL ROOTS 


Problems for Lecture 14 
1. Solve the following initial value problem: Y + 4% + 3x = 0, with x(0) = 1 and x(0) =0. 
2. Find the solution of Y — x = 0, with x(0) = xo and x(0) = uo. 


Solutions to the Problems 


Lecture 15 


Case 2: complex-conjugate roots (Part 
A) 


View this lecture on YouTube 


When the roots of the characteristic equation are complex conjugates, we can define real numbers 


A and u such that the two roots are given by 
r=A+ip, F=A-in. 
We have thus found the following two complex exponential solutions to the differential equation: 
z(t) = eit, ee, 


Applying the principle of superposition, any linear combination of z and Z is also a solution to the 
second-order ode, and we can form two different linear combinations of z(t) and Z(t) that are real, 
namely x1 (t) = Re z(t) and x2(t) = Im z(t). We have 


x1(t) =e™ cos pt, x(t) = e™ sin ut. 


Having found these two real solutions, we can then apply the principle of superposition a second time 


to determine the general solution for x(t): 
x(t) = e™ (A cos ut + Bsin ut). 


The real part of the roots of the characteristic equation appear in the exponential function, the imagi- 


nary part appears in the cosine and sine. 
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LECTURE 15. CASE 2: COMPLEX-CONJUGATE ROOTS (PART A) 


Lecture 16 


Case 2: complex-conjugate roots (Part 
B) 


View this lecture on YouTube 


Example: Solve <+x+x=0withx(0) = 1 and x(0) = 0. 


The characteristic equation is r+r+1= 0, with roots 


1 v3 _ 1 .y3 
= 


t= Hi r= 1 
2° 2 


The general solution of the ode is therefore 
x(t) =e" (a cos = + Bsin 2) : 
The derivative is 


x(t) = ze (Acos Ya + Bsin 2) ep Asn i el 2) 


Applying the initial conditions x(0) = 1 and x(0) = 0 results in 


with solution A = 1 and B = \/3/3. Therefore, the resulting solution is 


x(t) = e? (s G, H V ie 2) : 


3 2 
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48 LECTURE 16. CASE 2: COMPLEX-CONJUGATE ROOTS (PART B) 


Problems for Lecture 16 
1. Solve £—-2x+5x = 0, with x(0) = 1 and x(0) = 0. 
2. Solve &+x = 0, with x(0) = xy and x(0) = uo. 


Solutions to the Problems 


Lecture 17 


Case 3: Repeated roots (Part A) 


View this lecture on YouTube 


If the exponential ansatz yields only one linearly independent solution, we need to find the miss- 


ing second solution. We will do this by starting with the roots of the characteristic equation given by 


r= À + ip, and taking the limit as y goes to zero. 
Now, the general solution for the case of complex-conjugate roots is given by 


x(t) = e™ (Acos ut + Bsin ut), 


and to limit this solution as u — 0 requires first satisfying the general initial conditions x(0) = xy and 


x(0) = uo. Solving for A and B, the solution becomes 


Xo 


x(t; u) = e™ (x cos ut + ae sin ut) j 


where we have written x = x(t; u) to show explicitly that x now depends on the parameter u. 
Taking the limit as 4 > 0, and using lim, +o A sin ut = t, we have 
lim x(t; u) = e™ (xo + (uo — Axo)t); 


u—>0 


and the second solution is observed to be a constant, uy — Axo, times t times the first solution. The 
general solution for the case of repeated roots can therefore be written in the form 


x(t) = (c1 + ct)e", 


where r is the repeated root. The main result is that the second solution is t times the first solution. 


49 


50 


LECTURE 17. CASE 3: REPEATED ROOTS (PART A) 


Lecture 18 


Case 3: Repeated roots (Part B) 


View this lecture on YouTube 
Example: Solve <+2x+x=0withx(0) = 1 and x(0) = 0. 


The characteristic equation is 1? + 2r +1 = (r + 1)? = 0, which has a repeated root given by r = —1. 


Therefore, the general solution to the ode is 
x(t) = (cy +c2t)e*, x(t) = (co — c1 — cote. 


Applying the initial conditions, we have cj = 1 and cz — c1 = 0, so that cı = c2 = 1. Therefore, the 


solution is 
x(t) = (1 + t)e™. 
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52 LECTURE 18. CASE 3: REPEATED ROOTS (PART B) 


Problems for Lecture 18 
1. Solve £-2x+x=0,withx(0) = 1 and x(0) = 0. 


Solutions to the Problems 


Practice quiz: Homogeneous equations 
1. Solve & — 3X + 2x = 0, with initial values x(0) = 1 and x(0) = 0. 

2. Solve ¢ — 2% + 2x = 0, with initial values x(0) = 1 and x(0) = 0. 

3. Solve Y +2 + x = 0, with initial values x(0) = 0 and x(0) = 1. 


Solutions to the Practice quiz 
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LECTURE 18. CASE 3: REPEATED ROOTS (PART B) 


Lecture 19 


Inhomogeneous second-order ode 


View this lecture on YouTube 


The inhomogeneous linear second-order ode is given by 
ë+ p(t)x + q(t)x = g(t), 


with initial conditions x(to) = xy and x(tp) = uo. There is a three-step solution method when the 


inhomogeneous term g(t) # 0. (1) Solve the homogeneous ode 
ï+ p(t)x+q(t)x =0, 


for two independent solutions x = x;(t) and x = x(t), and form a linear superposition to obtain the 
general solution 
xp(t) = c1x1 (t) + cox2(t), 


where cı and cp are free constants. (2) Find a particular solution x = xp(t) that solves the inhomoge- 
neous ode. A particular solution is most easily found when p(t) and q(t) are constants, and when g(t) 
is a combination of polynomials, exponentials, sines and cosines. (3) Write the general solution of the 
inhomogeneous ode as the sum of the homogeneous and particular solutions, 


x(t) = xp(t) + xpt), 


and apply the initial conditions to determine the constants cı and cz appearing in the homogeneous 


solution. Note that because of the linearity of the differential equations, 


x d? d 
X+ pi + qx = 75 (xp + Xp) 4 Pa (Xn Xp) + (Xp + Xp) 


= (Ën + pXp +qXp) + (Xp + pXp + Xp) 
=0+8 
= 8, 


so that the sum of the homogeneous and particular solutions solve the ode, and the two free constants 
in x, can be used to satisfy the two initial conditions. 
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56 LECTURE 19. INHOMOGENEOUS SECOND-ORDER ODE 


Problems for Lecture 19 
1. Consider the inhomogeneous linear second-order ode given by 
¥+ p(t) +q(t)x = gi(t) + galt). 
Show that 
x(t) = xn (t) + Xp, (E) + Xp, (t) 


is the general solution, where x;,(t) is the general solution to the homogeneous ode, xp, (t) is a par- 
ticular solution for the inhomogeneous ode with only g(t) on the right-hand-side, and xp,(t) is a 
particular solution for the inhomogeneous ode with only g2(t) on the right-hand side. 


Solutions to the Problems 


Lecture 20 


Inhomogeneous term: exponential 


function 


View this lecture on YouTube 
Example: Solve & — 3x — 4x = 3e” with x(0) = 1 and x(0) = 0. 
First, we solve the homogeneous equation. The characteristic equation is 


r? —3r—4= (r—4)(r+1)=0, 


so that 


xp(t) = cie” tet, 


Second, we find a particular solution of the inhomogeneous equation. We try an ansatz such that the 


exponential will cancel out of both sides of the ode: 
x(t) = Ae”, 
where A is an undetermined coefficient. Substituting our ansatz into the ode, we obtain 
4A-6A-4A=3, 


from which we determine A = —1/2. Obtaining a solution for A independent of t justifies the ansatz. 
Third, we write the general solution to the ode as the sum of the homogeneous and particular solu- 
tions, and determine cı and cp that satisfy the initial conditions. We have 


1 
x(t) = cye™ +C — 2° , x(t) = 4c1e“ = Ce! =e, 


Applying the initial conditions yields the two linear equations 


ata=5, 441 —c2 = 1, 


with solution cı = 1/2 and cp = 1. Therefore, the solution for x(t) that satisfies both the ode and the 


initial conditions is given by 


1 1 1 
x(t) = ze ze Lats =e (1 et, =) ; 


where we have grouped the terms in the solution to better display the asymptotic behavior for large t. 
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58 LECTURE 20. INHOMOGENEOUS TERM: EXPONENTIAL FUNCTION 


Problems for Lecture 20 


1. Solve Y +5% + 6x = e *, with x(0) = 0 and x(0) = 0. 


Solutions to the Problems 


Lecture 21 


Inhomogeneous term: sine or cosine 
(Part A) 


View this lecture on YouTube 


Example: Find a particular solution of Y — 3x — 4x = 2sint. 
We will demonstrate two solution methods in this and the next lecture. The first method tries the 


ansatz 
x(t) = Acost+ Bsint, 


where the argument of cosine and sine must agree with the argument of sine in the inhomogeneous 
term. The cosine term is required because the derivative of sine is cosine. Upon substitution into the 


differential equation, we obtain 
(-Acost — Bsint) — 3 (—Asint + Bcost) —4(Acost+ Bsint) = 2sint, 


or regrouping terms, 
— (5A + 3B) cost + (3A — 5B) sin t = 2sint. 


This equation is valid for all t, and in particular for t = 0 and 7/2, for which the sine and cosine 


functions vanish. For these two values of t, we find 
5A+3B=0, 3A-5B=2; 
and solving we obtain A = 3/17 and B = —5/17. The particular solution is therefore given by 


1 
Xp = 75 (3cost — 5sint). 
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LECTURE 21. INHOMOGENEOUS TERM: SINE OR COSINE (PART A) 


Lecture 22 


Inhomogeneous term: sine or cosine 
(Part B) 


View this lecture on YouTube 

Example: Find a particular solution of Y — 3x — 4x = 2sint. 

The second method makes use of the relation e# = cost + isin t to convert the sine inhomogeneous 
term to an exponential function. We introduce the complex function z = z(t) and note that sint = 


Im(e**). We then consider the complex ode 
2 — 32 —4z = 2e", 


where x = Im{z} satisfies the original differential equation for x. 
To find a particular solution of the complex equation, we try the ansatz z(t) = Ce’, where we now 
expect C to be a complex constant. Upon substitution into the complex ode, and using i? = —1, we 


obtain —C — 3iC — 4C = 2; or upon solving for C: 


-2 543i 


C= +31 17 


Therefore, 
1 ; IR 1 l 
Xp = Im{zp} = Im { y (5 +31) (cost-+isins) | = 7 cost —5sint), 


the same result as in the first method. 
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62 LECTURE 22. INHOMOGENEOUS TERM: SINE OR COSINE (PART B) 


Problems for Lecture 22 


1. Find a particular solution of Y — 3x — 4x = 2cost by trying a 
a) cosine and sine ansatz; 
b) an exponential ansatz. 


Solutions to the Problems 


Lecture 23 


Inhomogeneous term: polynomials 


View this lecture on YouTube 
Example: Find a particular solution of &+x— 2x = t?. 
Our ansatz should be a polynomial in t of the same order as the inhomogeneous term. Accordingly, 


we try 
x(t) = At? + Bt +C. 


Upon substitution into the ode, we have 


2A + (2At+ B) —2(AÊ +Bt+C) = f, 


or 
2AP +2(A—B)t+ (2A +B -2C =P. 


Equating powers of t, 
—2A=1, 2(A-B)=0, 24+B-2C=0; 


and solving, 


The particular solution is therefore 
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64 LECTURE 23. INHOMOGENEOUS TERM: POLYNOMIALS 


Problems for Lecture 23 


1. Find a particular solution of ¥+%+x=t. 


Solutions to the Problems 


Lecture 24 


Resonance 


View this lecture on YouTube 
Resonance occurs when the forcing frequency matches the natural frequency of an oscillator. To 


illustrate resonance, we consider the inhomogeneous linear second-order ode 
¥+w hx =fcoswt, x(0) =0, x(0) =0; 


and determine what happens to the solution in the limit w — wo. 


The homogeneous equation has characteristic equation r2 + ws = 0, with solution r = tiwg, so 


that the solution to the homogeneous equation is 
xp(t) = c1 cos Wot + cz sin wot. 


A particular solution of the equation can be found by trying x(t) = A cos wt. Upon substitution into 
the ode, we can solve for A to obtain A = f/(w} — w°), so that 


xp(t) = zz a s wt. 
Our general solution is thus 


x(t) = c1 cos wot + c2 sin wot + a cos wt. 
we — w 


0 
Initial conditions are satisfied when cı = —f /(w4 — w?) and c = 0, so that 


cos wt — cos wot) 


x(t) z 


2.42 
ww 


Resonance occurs in the limit  — wo; that is, the frequency of the inhomogeneous term (the ex- 
ternal force) matches the frequency of the homogeneous solution (the free oscillation). By L'Hospital's 
rule, the limit of x = x(t) is found by differentiating with respect to w: 


—ftsinwt _ ftsinwot 


cos wt — cos wot 
lim x(t) = lim A OE fie 
WwW WwW ws — q? ww —2w 2wo 


At resonance, the oscillation increases linearly with t, resulting in a large amplitude oscillation. 
More generally, if the inhomogeneous term in the differential equation is a solution of the cor- 
responding homogeneous differential equation, then one should multiply the usual ansatz for the 


particular solution by t. 
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66 LECTURE 24. RESONANCE 


Problems for Lecture 24 


1. Solve # + 3x +2x = e~*!, with x(0) = 0 and x(0) = 0. Observe that the inhomogeneous term is a 
solution of the homogeneous equation. To find a particular solution, the usual exponential ansatz will 


need to be multiplied by t. 


Solutions to the Problems 


Practice quiz: Inhomogeneous 


equations 


1. Find a particular solution of Y + 3x + 2x = 202. 
2. Find a particular solution of Y — x — 2x = 2cos2f. 
3. Find a particular solution of Y —3x+2x = t + 1. 


Solutions to the Practice quiz 
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LECTURE 24. RESONANCE 


Lecture 25 


Application: RLC circuit 


View this lecture on YouTube 


e(t) 


RLC circuit diagram 


Consider a resister R, an inductor L, and a capacitor C connected in series as shown in the above figure. 
An AC generator provides a time-varying electromotive force to the circuit, given by E(t) = Eg cos wt. 
The equations for the voltage drops across a capacitor, a resister and an inductor are 
di 
Ve=g/C, Vr=iR, V= q 

where C is the capacitance, R is the resistance, and L is the inductance. The charge q and the current i 
are related by i = dq /dt. 

Kirchhoff's voltage law states that the electromagnetic force applied to any closed loop is equal to 
the sum of the voltage drops in that loop. Appying Kirchhoff's law, we have V + Vr + Ve = €(t), or 


2 
aq. p4q 1 


Pop et cd 


= Ey cos wt. (25.1) 


This is a second-order linear inhomogeneous differential equation with constant coefficients. 

To reduce the number of free parameters in this equation, we can nondimensionalize. We first 
define the natural frequency of oscillation of a system to be the frequency of oscillation in the absence 
of any driving or damping forces. For the RLC circuit, the natural frequency of oscillation is given 
by wo = 1/ VLC, and making use of wọ, we can define a dimensionless time T and a dimensionless 
charge Q by 


wel 
T=wt, Q= ar 
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70 LECTURE 25. APPLICATION: RLC CIRCUIT 


The resulting dimensionless equation for the RLC circuit can then be found to be 


PQ dQ 
da ae +Q = cos Br, (25.2) 


where « and ß are dimensionless parameters given by 


Lwo i wo ` 


Notice that the original five parameters in (25.1), namely R, L, C, Eq and w, have been reduced to the 
two dimensionless parameters « and ß in (25.2). 
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Problems for Lecture 25 


1. Determine how to nondimensionalize the RLC circuit equation so that the dimensionless equation 
takes the form 


PQ dQ, 
sa“ HQ = cos BT. 


What are the definitions of t, Q, x and 6? 


a 


Solutions to the Problems 


72 


LECTURE 25. APPLICATION: RLC CIRCUIT 


Lecture 26 


Application: mass on a spring 


View this lecture on YouTube 


l, m 
YA ES 


Mass-spring system 


Consider a mass connected by a spring to a fixed wall, as shown above. The spring force is modeled 

by Hooke's law, F; = —kx, and sliding friction is modeled as Ff = —cdx/dt. An external force is 

applied and is assumed to be sinusoidal, with F. = Fy cos wt. Newton’s equation, F = ma, results in 
dx dx 


ma ch Hkx = Focos wt, 


another second-order linear inhomogeneous differential equation with constant coefficients. Here, the 
natural frequency of oscillation is given by wo = vk/m, and we define a dimensionless time T and a 
dimensionless position X by 


w t X x 
T= ot, = A 
The resulting dimensionless equation is given by 
| a + X co pt ( 6 ) 
2 T S , 2 1 


where here, a and f are dimensionless parameters given by 


Even though the physical problem is different, the resulting dimensionless equation for the mass- 


spring system is the same as that for the RLC circuit. 
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74 LECTURE 26. APPLICATION: MASS ON A SPRING 


Problems for Lecture 26 


1. Determine how to nondimensionalize the mass on a spring equation so that the dimensionless 


equation takes the form 


What are the definitions of t, X, x and f? 


Solutions to the Problems 


Lecture 27 


Application: pendulum 


View this lecture on YouTube 


The pendulum 


Consider a mass attached to a massless rigid rod and constrained to move along an arc of a circle 
centered at the pivot point (see figure above). Suppose l is the fixed length of the connecting rod, and 
0 is the angle it makes with the vertical. 

We apply Newton's equation, F = ma, to the mass with origin at the bottom and axis along the 
arc with positive direction to the right. The position s of the mass along the arc is given by s = 10. 
The relevant gravitational force on the pendulum is the component along the arc, and is given by 
Fs = —mgsin@. We model friction to be proportional to the velocity of the pendulum along the 
arc, that is Fr = —c§ = —cl§. With a sinusoidal external force, F, = Fo cos wt, Newton’s equation 
ms = Fe + Fr + F; results in 

mlé + cl0 + mg sinó = Focos wt. 


At small amplitudes of oscillation, we can approximate sin ~ 0, and the natural frequency of 
oscillation wo of the mass is given by wo = y2/1. Nondimensionalizing time as tT = wot, the dimen- 
sionless pendulum equation becomes 

d?9 dé 


-57 tar + sin = y cos PT, 


where «, ß, and y are dimensionless parameters given by 


76 LECTURE 27. APPLICATION: PENDULUM 


The nonlinearity of the pendulum equation, with the term sin 0, results in the additional dimensionless 
parameter y. For small amplitude of oscillation, however, we can scale 0 by 0 = yO, and the small 
amplitude dimensionless equation becomes 

Po de 


T2 | a + © = cos BT, (27.1) 


the same equation as the dimensionless equations for the RLC circuit and the mass on a spring. 


Lecture 28 


Damped resonance 


View this lecture on YouTube 
The dimensionless odes for the RLC circuit, mass on a spring, and small amplitude pendulum have 
the form 

¥+ax+x=cos Bt, (28.1) 


where the physical constraints of these applications require that a > 0. The corresponding homoge- 
neous equation has the characteristic polynomial r? + ar +1 = 0, with roots r+ = (—a + va? —4)/2. 
No matter the sign of the discriminant, we have Re{r+} < 0. Therefore, both linearly independent 


homogeneous solutions decay exponentially to zero, and the long-time solution of the differential 
equation reduces to the non-decaying particular solution. Since the initial conditions are satisfied by 
the free constants multiplying the decaying homogeneous solutions, the long-time solution is inde- 
pendent of the initial conditions. 

If we are only interested in the long-time solution of (28.1), we only need determine the particular 


solution. We can consider the complex ode 
Z+az+z=clbt, 
with x, = Re(zp). With the ansatz zp = Ae'P!, we have —f?A + inBA + A = 1, or 


1 
(=) + inp" 


Of particular interest is when the forcing frequency w is equal to the natural frequency wo of the 


A= 


undampled oscillator. In the dimensionless equation of (28.1), this resonance situation corresponds to 
b = 1. In this case, A = 1/ix and 


Xp = Re{e"/ia} = (1/«) sint. 


The oscillator position is observed to be 71 /2 out of phase with the external force, or in other words, the 
velocity of the oscillator, not the position, is in phase with the force. Also, the amplitude of oscillation 
is seen to be inversely proportional to the damping coefficient represented by a. Smaller and smaller 


damping will obviously lead to larger and larger oscillations. 
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78 LECTURE 28. DAMPED RESONANCE 


Problems for Lecture 28 


1. Consider the differential equation 
X+ ax +x = cos Bt. 


Find the long-time amplitude of oscillation as a function of « and $£. 


Solutions to the Problems 


Week III 


The Laplace Iransform and Series 
Solution Methods 


81 


We present here two new analytical solution methods for solving linear differential equations. 
The first is the Laplace transform method, which is used to solve the constant-coefficient ode with a 
discontinuous or impulsive inhomogeneous term. The Laplace transform is a good vehicle in general 
for introducing sophisticated integral transform techniques within an easily understandable context. 
We also introduce the solution of a linear differential equation by series solution. Although we do not 
go deeply here, an introduction to this technique may be useful to students that encounter it again in 


more advanced courses. 
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Lecture 29 


Definition of the Laplace transform 


View this lecture on YouTube 
The Laplace transform of f(t), denoted by F(s) = L{f(t)}, is defined by the integral 
Els) = f tdt. (29.1) 


The values of s may be restricted to ensure convergence. The Laplace transform can be shown to be a 


linear transformation. We have 


apio Hefe} = f P e (c1 filt) + cofa(t)) at 
Cy I ef (tdt + cz [ ed 


=CL1fil(t)) +coL{fo(t)}. 


There is also a one-to-one correspondence between functions and their Laplace transforms, and a 
table of Laplace transforms is used to find both Laplace and inverse Laplace transforms of commonly 
occurring functions. 


To construct such a table, integrals have been performed, where for example, we can compute 


eier} = f° eas 
BE 1 ean 
s—a 0 
o 1 
ea 


Here, the assumption is made that s > a. 
A Table of Laplace transforms can be found in Appendix F. 
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84 LECTURE 29. DEFINITION OF THE LAPLACE TRANSFORM 


Problems for Lecture 29 


1. Compute the Laplace transform of f(t) = sin bt. 


Solutions to the Problems 


Lecture 30 


Laplace transform of a constant 


coefficient ode 


View this lecture on YouTube 


Consider the inhomogeneous constant-coefficient second-order ode, given by 
aš be Pee S g(t), x(0)= xo, ž(0) = uo. 
We Laplace transform the ode using linearity of the transform: 
al{&} +bL1x) +c£1x) = L1g). 


We define X(s) = L{x(t)} and G(s) = L{g(t)}. The Laplace transforms of the derivatives can be 
found by an integration by parts. We have 


oo 


oo : y sad c 
/ etxdt=xe | + s/ e!xdt = sX(s) — xo; 
0 0 


oo 


/ eidt=x0e | + s/ e“ żdt = —ug +s(sX(s) — xo) = s?X(s) — sxo — uo. 
0 0 0 


The transformed ode is then 


a (s?x SXxo uo) Hb (sX — xo) +cX =G, 


which is an easily solved linear algebraic equation for X = X (s). To finish the solution of the ode then 
requires taking the inverse Laplace transform of X = X (s) to obtain x = x(t). 


85 


86 LECTURE 30. LAPLACE TRANSFORM OF A CONSTANT COEFFICIENT ODE 


Problems for Lecture 30 
1. Consider the inhomogeneous constant-coefficient second-order ode given by 
ax Hoe per =e), x(0)= xo, ž(0) = uo. 


Determine the solution for X = X(s) in terms of the Laplace transform of g(t). 


Solutions to the Problems 


Lecture 31 


Solution of an initial value problem 


View this lecture on YouTube 
Example: Solve * + x = sin2t with x(0) = 2 and x(0) = 1 by Laplace transform methods. 
Taking the Laplace transform of both sides of the ode, we find 


2 

2 

X(s)—2s—1+X(s) = =—, 

EX) 25-14 X(5) = z 

where the Laplace transform of the second derivative and of sin 2t made use of line 16 and line 6 of 


the table in Appendix F. Solving for X(s), we obtain 


Gel. 2 
s?+1 (s2+1)(s2 +4) 


X(s) 


To determine the inverse Laplace transform from the table in Appendix F, we need to perform a partial 


fraction expansion of the second term: 


2 _as+b cs+d 


: 1.1 
(s2+1)(s?+4) s?+1 5244 oe) 


By inspection, we can observe that a = c = 0 and that d = —b. A quick calculation shows that 3b = 2, 
or b = 2/3. Therefore, 


Sri 2/8. 2/3 
2+1 32+1  (s2+4) 

© 2s 5/3 2/3 
2+1 2+1 (2+4) 


X(s) 


From lines 6 and 7 of the table in Appendix F, we obtain the solution by taking inverse Laplace 
transforms of the three terms separately, where the values in the table are b = 1 in the first two terms, 
and b = 2 in the third term: 

x(t) = 2cos t + sint — sin 2t. 
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88 LECTURE 31. SOLUTION OF AN INITIAL VALUE PROBLEM 


Problems for Lecture 31 


1. Solve Y +5% + 6x = e', with x(0) = 0 and x(0) = 0, by Laplace transform methods. 


Solutions to the Problems 


Practice quiz: The Laplace transform 
method 


1. Find the Laplace transform of x(t) = e”' cos rt. 
2. For + Xx — 6x =e"*, with x(0) = x(0) = 0, find X(s). 


1 
(s+1)(s+2)(s+3) 


3. For X(s) = ‚find x(t). 


Solutions to the Practice quiz 
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LECTURE 31. SOLUTION OF AN INITIAL VALUE PROBLEM 


Lecture 32 


The Heaviside step function 


View this lecture on YouTube 


The Heaviside or unit step function, denoted here by u.(t), is zero for t < c and one fort > c: 


0, t<c; 
u-(t) = 32.1 
e(t) e t>c. pad 


The Heaviside function can be viewed as the step-up function. Using the Heaviside function both 
a step-down and a step-up, step-down function can also be defined. The Laplace transform of the 


Heaviside function is determined by integration: 


Lay = Puga = 


The Heaviside function can be used to represent a translation of a function f(t) a distance c in the 
positive t direction. We have 

0, t<c; 
f(t-c), t>c. 


ult)f(t-c) = 


The Laplace transform is 


L{uc(t)f(t—e)} = I eu, (t) f(t - cjdi = e-F(s). 


The translation of f (t) a distance c in the positive t direction corresponds to the multiplication of F(s) 


= CS. 


by the exponential e 
Piecewise-defined inhomogeneous terms can be modeled using Heaviside functions. For example, 
consider the general case of a piecewise function defined on two intervals: 


_ fit), ift<c; 
an pl), ift>c. 


Using the Heaviside function uç, the function f (t) can be written in a single line as 


f(t) = ft) + (2) — fi) u(t). 
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92 LECTURE 32. THE HEAVISIDE STEP FUNCTION 


Problems for Lecture 32 
1. Using the step-down function u.(t), construct a step-up and a step-up, step-down function. 
2. Prove that L{u.(t)f(t-c)} =e “SF(s). 
3. Consider the piecewise continuous function given by 
fü) = k A < 1; 
1, ift>1. 
a) Express f(t) in a single line using the Heaviside function. 
b) Find F = F(s). 


Solutions to the Problems 


Lecture 33 


The Dirac delta function 


View this lecture on YouTube 


The Dirac delta function, denoted as ö(t), is defined by requiring that for any function f(t), 


[rwswat = fo. 


The usual view of the Dirac delta function is that it is zero everywhere except at t = 0, at which it is 
infinite in such a way that the integral is one. The Dirac delta function is technically not a function, 
but is what mathematicians call a distribution. Nevertheless, in most cases of practical interest, it can 
be treated like a function, where physical results are obtained following a final integration. 

There are several common ways to represent the Dirac delta function as a limit of a well-defined 
function. For our purposes, the most useful representation of the shifted Dirac delta function, 5(f — c), 


makes use of the step-up, step-down function constructed from Heaviside functions: 


ö(t- c) = lim ltc-elt) -Uere). 


Before taking the limit, the well-defined step-up, step-down function is zero except over a small inter- 
val of width 2e centered at t = c, over which it takes the large value 1/2e. The integral of this function 
is one, independent of the value of e. 

The Laplace transform of the Dirac delta function is easily found by integration using the definition 
of the delta function. With c > 0, 


ehe) I eE —c)dt = e7“, 
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94 LECTURE 33. THE DIRAC DELTA FUNCTION 


Problems for Lecture 33 
Remember that the Dirac delta function is well-defined only after it is integrated over. 


1. Prove that ö(ax) = ai) 


2. Show that 6,(x) = J d(x' —c)dx’. 


3. Show that d(x — c) = q de (*) : 


Solutions to the Problems 


Lecture 34 


Solution of a discontinuous 


inhomogeneous term 


View this lecture on YouTube 
Example: Solve % + 3x+2x =1—u1(t), with x(0) = x(0) = 0 


Here, the inhomogeneous term is a step-down function, from one to zero. Taking the Laplace trans- 


form of the ode using the table in Appendix F, we have 


s? X(s) +3sX(s) +2X(s) = 


ulm 


(=e), 
with solution for X = X (s) given by 


1-e”° 


N ner) 


(34.1) 


Defining 
1 
ET ENG 


and using the table in Appendix F, the inverse Laplace transform of X (s) can be written as 


x(t) = f(t) -—m (Hf E- 1), 


where f(t) is the inverse Laplace transform of F(s). To determine f(t), we need a partial fraction 


expansion of F(s), so we write 


1 Ooa, b > c 
s+2 


s(s+1)(s+2) s s+1 


Using the coverup method, we find a = 1/2, b = —1, and c = 1/2, and from the table in Appendix F, 


determine 


96 LECTURE 34. SOLUTION OF A DISCONTINUOUS INHOMOGENEOUS TERM 


Problems for Lecture 34 


1. Show that the solution in the lecture, 


is continuous at t = 1. 
2. Solve&+x=1- uzz (t), with x(0) = 0 and x(0) = 0. 


Solutions to the Problems 


Lecture 35 


Solution of an impulsive 


inhomogeneous term 


View this lecture on YouTube 
Example: Solve & + 3x + 2x = ö(t) with x(0) = x(0) = 0. Assume the entire impulse occurs at t = 0*. 


Taking the Laplace transform of the ode using the table in Appendix F, and applying the initial 


conditions, we have 
s*X + 3sK +2X =1. 


Solving for X = X(s), we have 
1 


%6) = ne 


and a partial fraction expansion results in 


1 1 1 


(s+1)(s+2) s+1 s+2 


Taking the inverse Laplace transform using the table in Appendix F, we find 


== 


=e"'(1-e"). 


Note that the function x = x(t) is continuous at t = 0, but that the first derivative is not, since 


x(07) = 0 and x(0*) = 1. Impulsive forces result in a discontinuity in the velocity. 
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98 LECTURE 35. SOLUTION OF AN IMPULSIVE INHOMOGENEOUS TERM 


Problems for Lecture 35 
1. Solve # +x = (t) — ö(t — 272), with x(0) = 0 and x(0) = 0. 


Solutions to the Problems 


Practice quiz: Discontinuous and 


impulsive inhomogeneous terms 


1. Define the above pictured function using Heaviside step functions. 
2. Find the solution to # + x = 1 — uzr (t), with x(0) = 1 and x(0) = 0. 
3. Find the solution to &+x = 6(t) — 9(t— 2rr), with x(0) = 1 and x(0) = 0. 


Solutions to the Practice quiz 
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100 LECTURE 35. SOLUTION OF AN IMPULSIVE INHOMOGENEOUS TERM 


Lecture 36 


The series solution method 


View this lecture on YouTube 


Example: Find the general solution of y” + y = 0. 


View tutorial on YouTube 


By now, you should know that the general solution is y(x) = agcosx + a1 sin x, with ag and a1 con- 


stants. To find a power series solution, we write 
foe) 
y(x) = y anx"; 
n=0 
and upon differentiating term-by-term 
00 oo 
y(x) = $} nae, Ve) nad 
n=1 
Substituting the power series for y and its derivatives into the differential equation, we obtain 
foe} oo 
Y n(n—1)anx"?+ } anx" =0. (36.1) 
n=2 n=0 


The power-series solution method requires combining the two sums into a single power series in x. 


We shift the summation index downward by two in the first sum to obtain 


E n(n Dana = Tin +2) (n + 1)an+2X”. 


We can then combine the two sums in (36.1) to obtain 


oo 


L (o F 2)(n Fe Dan+2 + an)" =0. 
n=0 
For the equality to hold, the coefficient of each power of x must vanish separately. We therefore obtain 


the recurrence relation 
An 


an2 = eat 


We observe that even and odd coefficients decouple. We thus obtain two independent sequences 


n=0,1,2,.... 
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102 LECTURE 36. THE SERIES SOLUTION METHOD 


starting with first term ag or a1. Developing these sequences, we have for the first sequence, 


1 1 
10, MT 44 = ~ 4342 g 
and so on; and for the second sequence, 
a a3 = — 1 a a5 = l az = —a 
1, 3S 3.2 1, 5 5.4 Be 1, 


and so on. Using the principle of superposition, the general solution is therefore 


Karls HBS | ob 
v(x) =a (1 T | AT =) En (x 31 | 51 =) 


= ag cos x + a, sin x, 


as expected. 
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Problems for Lecture 36 


1. Find two independent power series solutions of y” — y = 0. Show that 
y(x) = ay cosh x + ay sinh x, 


where 


e 
cosh x = 


Solutions to the Problems 


104 LECTURE 36. THE SERIES SOLUTION METHOD 


Lecture 37 


Series solution of the Airy’s equation 
(Part A) 


View this lecture on YouTube 


We solve by series solution the Airy’s equation, an ode that arises in optics, fluid mechanics, and 
quantum mechanics. 
Example: Find the general solution of y” — xy = 0. Notice that there is a non-constant coefficient. 


We try the power-series ansatz 
oo 
y(x) = Y anx", 
n=0 


where the a, coefficients are to be determined. Finding the second derivative by differentiating term- 


by-term, the ode becomes 
oo 


P n(n—1)anx"™? — Y anx"! =0. (37.1) 
n=2 n=0 


We shift the summation index of the first sum down by three to obtain 


00 


3 n(n—-1)anx"?= Y (n+3)(n + 2)an43x"*!. 


n=-1 


a 
ll 
N 


Y 


When combining the two sums in (37.1), we separate out the extra first term in the first sum. Therefore, 


(37.1) becomes the single power series 
209 + Y (o +3) (n +2)a,+3 — an) m0, 
n=0 


Since all the coefficients of powers of x equal zero, we find a2 = 0 and obtain the recursion relation 


1 
n+3)(n +2 


An+3 = ( jr 


Three sequences of coefficients—those starting with either ag, aı or aa—decouple. Since a2 = 0, we 


find immediately that a2 = a5 = ag = 411 = --- = 0. Starting with ag, we have 


1 1 
rm 


106 LECTURE 37. SERIES SOLUTION OF THE AIRY’S EQUATION (PART A) 


and so on; and starting with a4, 


= a a7 = l 
E E a 


and so on. The general solution for y = y(x), can therefore be written as 


1 x x6 x4 x7 
u(x) = 00 ( "Ba 6255852. =) a (x "2.377.643" de) 


= agyo(x) + a1y1 (x). 
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Problems for Lecture 37 


1. Find two independent power series solutions of y” + xy’ — y = 0. Keep terms up to x°. 


Solutions to the Problems 


108 LECTURE 37. SERIES SOLUTION OF THE AIRY’S EQUATION (PART A) 


Lecture 38 


Series solution of the Airy’s equation 
(Part B) 


View this lecture on YouTube 


The Airy’s equation is given by y” — xy = 0, and the general solution for y = y(x), is given by 


x3 x6 x4 x7 
vr) (tt ta (trat) 


= agyo(x) + a1y1 (x). 


Suppose we would like to graph the solutions y = yo(x) and y = yı(x) versus x by solving the 
differential equation y” — xy = 0 numerically. What initial conditions should we use? Evidently, 
y = yo(x) solves the ode with initial values y(0) = 1 and y’(0) = 0, whereas y = yı(x) solves the ode 
with initial values y(0) = 0 and y'(0) = 1. 

The numerical solutions, obtained using MATLAB, are shown in the figure below. Note that the 
solutions oscillate for negative x and grow exponentially for positive x. This can be understood by 
recalling that y” + y = 0 has oscillatory sine and cosine solutions and y” — y = 0 has exponential 


solutions. 


Airy’s functions 


Yo 
o 


Y, 
o 


-10 -8 -6 4 -2 0 2 


Numerical solution of Airy’s equation 
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110 LECTURE 38. SERIES SOLUTION OF THE AIRY’S EQUATION (PART B) 


Problems for Lecture 38 


1. Given the two Airy’s functions y = yo(x) and y = yı(x) as defined by their power series, solve the 
Airy’s equation y” — xy = 0 with the initial conditions y(0) = 1 and y'(0) = 1. 


Solutions to the Problems 


o e o o 
Practice quiz: Series solutions 
1. What does cosh? t — sinh? t equal? 

2. Find the general solution to y” + x?y = 0 up to terms proportional to x°. 
4 


3. Find the general solution to y” — xy’ + y = 0 up to terms proportional to x*. 
8 y yry P prop 


Solutions to the Practice quiz 
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112 LECTURE 38. SERIES SOLUTION OF THE AIRY’S EQUATION (PART B) 


Week IV 


Systems of Differential Equations and 


Partial Differential Equations 


113 
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We solve a coupled system of homogeneous linear first-order differential equations with constant 
coefficients. This system of odes can be written in matrix form, and we explain how to convert 
these equations into a standard matrix algebra eigenvalue problem. We then discuss the important 
application of coupled harmonic oscillators and the calculation of normal modes. The normal modes 
are those motions for which the individual masses that make up the system oscillate with the same 
frequency. Next, to prepare for a discussion of partial differential equations, we define the Fourier 
series of a function. Then we derive the well-known one-dimensional diffusion equation, which is 
a partial differential equation for the time-evolution of the concentration of a dye over one spatial 
dimension. We proceed to solve this equation for a dye diffusing length-wise within a finite pipe. 


116 


Lecture 39 


Systems of homogeneous linear 


first-order odes 


View this lecture on YouTube 


We consider a system of homogeneous linear odes with constant coefficients given by 
X1 =axı +bx2, K2=cxı +dx2, 


which can be written in matrix form as 


or more succinctly as x = Ax. 
We take as our ansatz x(t) = ve", where v and A are independent of t. Upon substitution into the 


ode, we obtain 
Ave = Ave% ; 


and upon cancellation of the exponential, we obtain the eigenvalue problem 
Av = Av. 
The characteristic equation for our two-by-two matrix is given by 


det (A — AI) = A? — (a+d)A + (ad — bc) = 0. 


We will demonstrate the solution for two separate cases: (i) eigenvalues of A are distinct and real; 
(ii) eigenvalues of A are complex conjugates. The third case for which the eigenvalues are repeated 
will be omitted. These three cases are analogous to those previously considered when solving the 


homogeneous linear constant-coefficient second-order ode. 
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118 LECTURE 39. SYSTEMS OF LINEAR FIRST-ORDER ODES 


Problems for Lecture 39 


1. Consider the system of homogeneous linear odes with constant coefficients given by 
xy = axı +CxX>, ž = CX, + dx. 


Prove that the eigenvalues of the resulting characteristic equation are real. 


Solutions to the Problems 


Lecture 40 


Distinct real eigenvalues 


View this lecture on YouTube 


Example: Find the general solution of žų =X1 + x2, Xp = 4x1 + x2. 


The equation in matrix form is given by 


s-h DE 


or X = Ax. With our ansatz x(t) = ve^, we obtain the eigenvalue problem Av = Av. The characteristic 


equation of A is given by 


det (A — AI) = 4? — 24 -3 = (A—3)(A+1) = 0, 


and the eigenvalues are found to be A; = —1 and Az = 3. To determine the corresponding eigenvec- 


tors, we substitute the eigenvalues successively into (A — Al)v = 0, and denote the eigenvectors as 
T 


v= (on 0) and v? = (vn v22) . 
For A; = —1, and unknown eigenvector v1, we obtain from (A — Aıl)vı = 0 the result v2] = —2711. 
Recall that an eigenvector is only unique up to multiplication by a constant: we may therefore take 
vı = 1. For Az = 3, and unknown eigenvector v2, we obtain from (A — AzI)v2 = 0 the result 
v22 = 2017. Here, we take 017 = 1. 
Therefore, our eigenvalues and eigenvectors are given by 


Aı=-l, vi = : ER RE, 
Sy NAT , ERN E x 
1 1 2 2 2 2 


Using the principle of superposition, the general solution to the ode is given by 
x(t) = cyvyer! + cove! 


or explicitly writing out the components, 


xilt) = ce ™t + cee, xa(t) = —2cye* + 2cpe*. 
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120 LECTURE 40. DISTINCT REAL EIGENVALUES 


Problems for Lecture 40 
1. Using matrix algebra methods, find the general solution of žų = —x2 and X = —2x1 — Xp. 


Solutions to the Problems 


Lecture 41 


Complex-conjugate eigenvalues 


View this lecture on YouTube i i 
Example: Find the general solution of x1 = 27% +x and ž = —X1 — =X. 


2 
d x1 = —5 1 X1 
dt X2 —1 -4 X2 l 


The ansatz x = ve^ leads to the characteristic equation 


The equations in matrix form are 


det(A AI) = +A+2 =0, 


which has complex-conjugate solutions. We denote the two eigenvalues as 


1 = 1 


The eigenvectors also occur as a complex-conjugate pair, and the eigenvector v associated with the 


T 
eigenvalue A satisfies —iv; + v2 = 0. Normalizing with vı = 1, we have v = (1 i) . We have 
therefore determined two independent complex solutions to the ode, that is, 


ve! and vet, 


and we can form a linear combination of these two complex solutions to construct two independent 


real solutions. Namely, the two real solutions are Refve**) and Im{ve?!}. We have 


Re{veX!} = Re P er — px > ) 
! — sint 
At IN 349: _1, [sint 
Im{ve"'} = Im ea Be l 
{ cost 


Taking a linear superposition of these two real solutions yields the general solution to the ode, given 
by 
Lg cost sin t 
x=e Y|A . +B A 
— sint cost 
or explicitly writing out the components, 


x= e72! (Acost+Bsint), X2 = e72! (—Asint + Bcost). 


121 


122 LECTURE 41. COMPLEX-CONJUGATE EIGENVALUES 


Problems for Lecture 41 
1. Using matrix algebra methods, find the general solution of žų = xı — 2x2 and X2 = x1 + x2. 


Solutions to the Problems 


Practice quiz: Systems of differential 


equations 


1. Find the second-order ode that has the same characteristic equation as the system of odes given by 


žı = axı + bx and Xo = cx, + dan. 
2. Find the general solution of X = xı + 2x2 and x2 = 2x1 + X2. 
3. Find the general solution of žų = —2x1 + x2 and %2 = —xı — 2x2. 


Solutions to the Practice quiz 
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124 LECTURE 41. COMPLEX-CONJUGATE EIGENVALUES 


Lecture 42 


Coupled oscillators 


View this lecture on YouTube 


> | 5,» 
X} X, 
Coupled harmonic oscillators: two masses, three springs 


The normal modes of a physical system are oscillations of the entire system that take place at a single 
frequency. Perhaps the simplest example of a system containing two distinct normal modes is the 
coupled mass-spring system shown above. We will see that an eigenvector analysis of this system, 
similar to what we have just done for coupled first-order odes, will exhibit the normal modes and 
reveal the true nature of the oscillation. 

In the figure, the position variables x; and x2 are measured from the equilibrium positions of the 
masses. Hooke's law states that the spring force is linearly proportional to the extension length of the 
spring, measured from equilibrium. By considering the extension of the spring and the sign of the 


force, we write Newton’s law F = ma separately for each mass: 


mı = —kxı — K(x1 — x2), 


mx) = —kx2 — K(x2 — x1). 
Further rewriting by collecting terms proportional to xı and x2 yields 


mx, = —(k+ K)x, + Kxo, 
mx) = Kx, — (k + K)x2. 


The equations for the coupled mass-spring system form a system of two homogeneous linear second- 


order odes. In matrix form, mx = Ax, or explicitly, 
= d? xa Y [- (k+K) K x 
dl Xx) K (k+K)) A x2) 
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126 LECTURE 42. COUPLED OSCILLATORS 


Problems for Lecture 42 


/ Kk k 
A 


VA X1 X2 


Coupled harmonic oscillators: two masses, two springs 


1. Consider the mass-spring system shown above. Determine the matrix equation represented as 
mx = Ax. 


Solutions to the Problems 


Lecture 43 


Normal modes (eigenvalues) 


View this lecture on YouTube 
We now solve the two masses, three springs system with governing equation 


d2 X1 _ —(k + K) K X1 
META x2 u K —(k+K) x2 ` 


We try the ansatz x = ve't, and obtain the eigenvalue problem Av = Av, with A = mr?. The eigenvalues 


are determined by solving the characteristic equation 
det (A — AI) = (A +k + K} — K? =0; 
and the two solutions for A are given by 


Ay=—k, Ap = —(k +2K). 


The corresponding values of r in our ansatz x = vet, with r = £\V/A/m, are 
ry =ivk/m, MA, Tr = iy/(k+2K)/m, Ta. 


Since the values of r are pure imaginary, we know that xı(t) and x>(t) will oscillate with angular 
frequencies w1 = Im{rı} and wz = Im{ra}, that is, 


wı = vk/m, w= ı/(k+2K)/m. 


These are the so-called frequencies of the two normal modes. In the next lecture, we find the eigen- 


vectors and physically intepret the results. 
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128 LECTURE 43. NORMAL MODES (EIGENVALUES) 


Lecture 44 


Normal modes (eigenvectors) 


View this lecture on YouTube 
We are solving the two masses, three springs system with governing equation 


d2 X1 o —(k+K) K X1 
META x2 u K —(k+K) x2 ` 


We tried the ansatz x = ve, and obtained the eigenvalue problem Av = Av, with A = mr?. The 


eigenvalues were found to be 
Aı=-k, A= —(k+2K). 


The eigenvectors, or so-called normal modes of oscillations, correspond to the mass-spring system 


oscillating at a single frequency. The eigenvector with eigenvalue A, satisfies 
—Kv11 + Kv12 = 0, 


so that v11 = v12. The normal mode with frequency wı = vk/m thus follows a motion where x1 = x2. 
During this motion the center spring length does not change, and the frequency is independent of K. 
Next, we determine the eigenvector with eigenvalue Az: 


Kv] + Kim = 0, 


so that v,1 = —vV22. The normal mode with frequency un = y (k +2K)/m thus follows a motion 
where xı = —x2. During this motion the two equal masses symmetrically push or pull against each 
side of the middle spring. 

A general solution for x(t) can be constructed from the eigenvalues and eigenvectors. Our ansatz 
was x = ve”, and for each of two eigenvectors v, we have a pair of complex conjugate values for r. 
Accordingly, we first apply the principle of superposition to obtain four real solutions, and then apply 
the principle again to obtain the general solution. With wi = Yk/m and w = y(k+2K)/m, the 
general solution is given by 


1 1 
(a) = G) (A cos wt + B sin w1t) + ( ) (Ccos wat + Dsinwyt), 
X2 = 


where the now real constants A, B, C, and D can be determined from the four independent initial 
conditions, x1 (0), x2(0), X1(0), and x2(0). 
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130 LECTURE 44. NORMAL MODES (EIGENVECTORS) 


Problems for Lecture 44 


Y k k 
Coupled harmonic oscillators: two masses, two springs 


1. Consider the matrix equation obtained from the mass-spring system shown above. Find the angular 
frequencies and the eigenvectors of the two normal modes. 


Solutions to the Problems 


Practice quiz: Normal modes 


1. Determine the matrix equation for the above mass-spring system. 
2. Determine the angular frequencies of the normal modes for the above mass-spring system. 


3. Determine the eigenvectors of the normal modes for the above mass-spring system. 


Solutions to the Practice quiz 
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132 LECTURE 44. NORMAL MODES (EIGENVECTORS) 


Lecture 45 


Fourier series 


View this lecture on YouTube 


Fourier series will be needed to solve the diffusion equation. A periodic function f(x) with period 2L 


can be represented as a Fourier series in the form 


(45.1) 


a = NT . 
f(x) = D+ Y (ancos a + by sin 


n=1 


Determination of the coefficients ag,a1,a2,... and b1,b2,b3,... makes use of orthogonality relations 


for sine and cosine. We first define the Kronecker delta dy, as 


Sop 1, fn=m; 
nm 0, otherwise. 


The orthogonality relations for n and m positive integers are then given as the integration formulas 


F cos (=) cos(“** )dx = Long; F sin (=) sin (== ax = Edin, 


[© (=) sin (“= )ax = 0. 


To determine the coefficient an, we multiply both sides of (45.1) by cos (nrıx/L), change the dummy 
summation variable from n to m, and integrate over x from —L to L. To determine the coefficient 


bn, we do the same except multiply by sin (nrıx/L). To illustrate the procedure, when finding by, 


integration will result in 


co L 
L bm | sin "T sin = dx 


m=1 


Fl scosin a dx = 


= y bm (Loum) 


m=1 


= Lb„, 


which can be solved for bn. The final results for the coefficients are given by 


1 yl NTCX 1 ft | ATX 
An = zf fe) cos ax, bn = TRO sın ax. 
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134 LECTURE 45. FOURIER SERIES 


Problems for Lecture 45 


1. Let 


a = NTTX . NTX 
fx) =F + Y (ancos 7 + by sin i 


a) Show that f(x + 2L) = f(x), that is, f(x) is a periodic function with period 2L. 
b) Show that ay is twice the average value of f(x). 


Solutions to the Problems 


Lecture 46 


Fourier sine and cosine series 


View this lecture on YouTube 


We now know that the Fourier series of a periodic function with period 2L is given by 


with 
il if, Fe cos =p da, m=z ff) 


The Fourier series simplifies if f(x) is an even function such that f(—x) = f(x), or an odd function 
such that f(-x) = —f(x). We make use of the following facts. The function cos (nrıx/L) is an even 
function and sin (n7tx/L) is an odd function. The product of two even functions is an even function. 
The product of two odd functions is an even function. The product of an even and an odd function is 


an odd function. And if g(x) is an even function, then 


[star =2 [sou 


a g(x)dx =0. 


—L 


and if g(x) is an odd function, then 


Regarding the Fourier series, then, if f(x) is even we have 


Ay a f(x ) cos “dx, Da = 0; 


and the Fourier series for an even function is given by the Fourier cosine series 


If f(x) is odd, then 
an =0, T f(x ) sin dx; 


and the Fourier series for an odd function is given by the Fourier sine series 


x) = } by sin Tr, f(x) odd. 
n=1 
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136 LECTURE 46. FOURIER SINE AND COSINE SERIES 


Problems for Lecture 46 


1. Show that the Fourier series for an odd function satisfies f(0) = 0 and that the Fourier series for an 


even function satisfies f’(0) = 0. 


Solutions to the Problems 


Lecture 47 


Fourier series: example 


View this lecture on YouTube 


Example: Determine the Fourier series of the triangle wave, shown in the following figure: 


IA N 
Te 


-27 -T 0 T on 


The triangle wave 


Evidently, the triangle wave is an even function of x with period 271, and its definition over half a 
period is 


=1-— 
f(x) a 0O<x<r 


Because f(x) is even, it can be represented by a Fourier cosine series (with L = 71) given by 
aca cos nx, with a == |" fía) cosnras 
2 = n , n 1 0 . 


The coefficient a9 /2 is the average value of f(x) over one period, which is clearly zero. The coefficients 


An for n > 0 are 


2 f7 2 4 8/(n?n?), ifn odd; 
a == | (1-2) cosneax = z „cum =] he ee WO 
m Jo TT nr 


0, if n even. 


The Fourier cosine series for the triangle wave is therefore given by 


8  COSI3x  cos5x | 
f(x) = = (cosx + tt): 


Convergence of this series is rapid. As an interesting aside, evaluation of this series at x = 0, using 


f (0) = 1, yields in an infinite series for 777, which is usually written in the form 


137 


138 LECTURE 47. FOURIER SERIES: EXAMPLE 


Problems for Lecture 47 


1. Determine the Fourier series of the square wave, shown in the following figure: 


-27 -T 0 T on 


The square wave 


Solutions to the Problems 


Practice quiz: Fourier series 


1. The Fourier series of a periodic function of period 2L is given by 


_mM y ATX . NTX : 
fix)= z +2 (an cos G— + by sin — F Determine the value of f(0). 


-1 
27 -r 0 7 27 


2. Determine the Fourier series of the square wave in the above figure. 


-1 
27 -r 0 7 27 


3. Determine the Fourier series of the square wave in the above figure. 


Solutions to the Practice quiz 
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Lecture 48 


The diffusion equation 


View this lecture on YouTube 


Derivation of the diffusion equation 


To derive the diffusion equation in one spacial dimension, we imagine a still liquid in a pipe of constant 
cross sectional area. A small quantity of dye is placed uniformly across a cross section of the pipe and 
allowed to diffuse up and down the pipe. We define u(x,t) to be the concentration (mass per unit 
length) of the dye at position x along the pipe, and our goal here is to find the partial differential 
equation satisfied by u. 

At time t, the mass of dye in an infinitesimal pipe volume located between x and x + Ax is given 
to order Ax by 

M = u(x,t)Ax. 


This mass of dye diffuses into and out of the cross sectional ends as shown in the figure. We assume 
the rate of diffusion is proportional to the concentration gradient, a relationship known as Fick’s law 
of diffusion. Fick’s law of diffusion assumes the mass flux J (mass per unit time) across a cross section 
of the pipe is given by 

J = —Dux, 


where D > 0 is the diffusion constant, and we have used the notation uy = du/dx. The mass flux 
is opposite in sign to the gradient of concentration so that the flux is from high concentration to low 
concentration. The time rate of change in the mass of dye is given by the difference between the 
mass flux into and the mass flux out of the infinitesimal cross sectional volume. A positive mass flux 
signifies diffusion from left to right. Therefore, the time rate of change of the dye mass is given by 


dM 


NEN — (x + Ax, t); 
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142 LECTURE 48. THE DIFFUSION EQUATION 
or rewriting in terms of u(x,t), 

ur(x,t)Ax = D (ux(x,t) — ux(x + Ax,t)). 
Dividing by Ax and taking the limit Ax — 0 results in the diffusion equation, given by 


uz = Duy. 
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Problems for Lecture 48 


1. Consider the one-dimensional diffusion equation in a pipe of length L. Nondimensionalize the 
diffusion equation using L as the unit of length and L?/D as the unit of time. 


Solutions to the Problems 


144 LECTURE 48. THE DIFFUSION EQUATION 


Lecture 49 


Solution of the diffusion equation 


(separation of variables) 


View this lecture on YouTube 


We begin to solve the diffusion equation 
ur = Duxy 


for the concentration u(x,t) in a pipe of length L. The solution method we use is called separation of 
variables. We assume that u(x,t) can be written as a product of two other functions, one dependent 


only on position x and the other dependent only on time t. That is, we make the ansatz 
u(x,t) = X(x)T(E). 


We will find all possible solutions of this type—there will in fact be an infinite number—and then use 
the principle of superposition to combine them to construct a general solution. 
Substitution of our ansatz into the diffusion equation results in 


XT = DX"T, 
which we rewrite by separating the x and t dependence to opposite sides of the equation: 


x 1 T’ 
x DT 
The left-hand side of this equation is independent of t and the right-hand side is independent of x. 
Both sides of this equation are therefore independent of both t and x and must be equal to a constant, 
which we call —A. Introduction of this separation constant results in the two ordinary differential 
equations 
X"+AX=0, T'+ADT=0. (49.1) 


To proceed further, boundary conditions need to be specified at the pipe ends. We will assume here 
that the ends open up into large reservoirs of clear fluid so that u(0,t) = u(L,t) = 0 for all times. 


Applying these boundary conditions to our ansatz results in 
(0,1) =X(0)T(£) =0, u(L,t) = X(L)T(t) =0. 


Since these boundary conditions are valid for all t, we must have X(0) = X(L) = 0. These are called 
homogeneous Dirichlet boundary conditions. In the next lecture, we proceed to solve for X = X(x). 
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146 LECTURE 49. SOLUTION OF THE DIFFUSION EQUATION (SEPARATION OF VARIABLES) 


Problems for Lecture 49 


1. Ifthe ends of the pipe were closed preventing a further diffusion of dye, then the mass flux through 
the pipe’s ends would be zero, and the appropriate boundary conditions would be u,(0,t) = ux(L,t) = 
0 for all times. Determine the appropriate boundary conditions on X = X(x). 


Solutions to the Problems 


Lecture 50 


Solution of the diffusion equation 


(eigenvalues) 


View this lecture on YouTube 


We now solve the equation for X = X(x) with homogeneous Dirichlet boundary conditions: 
X”+AX=0, X(0)=X(L)=0. 


Clearly, the trivial solution X(x) = 0 is a solution, and we will see that nontrivial solutions exist only 
for discrete values of A. These discrete values and the corresponding functions X = X(x) are called 
the eigenvalues and eigenfunctions of the ode. 

The form of the general solution of the ode depends on the sign of A. One can show that nontrivial 
solutions exist only when A > 0. We therefore write A = u?, and determine the general solution of 


X" +X =0 


to be 
X(x) = Acos px + Bsin px. 


The boundary condition at x = 0 results in A = 0, and the boundary condition at x = L results in 
BsinuL = 0. 
The solution B = 0 results in the unsought trivial solution. Therefore, we must have 
sin uL = 0, 


which is an equation for u, with solutions given by yn = nrt/L, where n is a nonzero integer. We have 


thus determined all the unique eigenvalues A = u? > 0 to be 
A,=(nn/L)%, n=1,2,3,..., 


with corresponding eigenfunctions 
Xn = sin (n7x/L). 


There is no need here to include the multiplication by an arbitrary constant, which will be added later. 
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148 LECTURE 50. SOLUTION OF THE DIFFUSION EQUATION (EIGENVALUES) 


Problems for Lecture 50 


1. Show that the equations 
X”+AX=0, X(0)=X(L)=0 


have no nontrivial solutions for A < 0. 


2. Solve the following equation for X = X(x) with given boundary conditions: 
X” +AX=0, X’(0) = X’(L) =0. 


Solutions to the Problems 


Practice quiz: Separable partial 


differential equations 


1. In the diffusion equation u; = Du,,, determine the units of the diffusion constant D. 


2. The one-dimensional wave equation is given by uş = c?uxy. With u(x,t) = X(x)T(t), determine 


the two ordinary differential equations after separation. 


3. Determine the eigenvalues and eigenfunctions of the differential equation X” + AX = 0, with mixed 
boundary conditions X(0) = 0 and X'(L) =0. 


Solutions to the Practice quiz 
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150 LECTURE 50. SOLUTION OF THE DIFFUSION EQUATION (EIGENVALUES) 


Lecture 51 


Solution of the diffusion equation 


(Fourier series) 


View this lecture on YouTube 


With the eigenvalues A, = (n7t/L)*, the differential equation for T = T(t) becomes 
T + (n272D/1?) T =0, 


which has solution proportional to 
T x ¿A Dt/ L 
n= i 


Therefore, with ansatz u(x,t) = X(x)T(t) and eigenvalues A, we conclude that the functions 


Un(x,t) = sin (nmx/L)e™ ™ DEL 
satisfy the diffusion equation and the spatial boundary conditions for every positive integer n. 

The principle of linear superposition for homogeneous linear differential equations then states that 
the general solution to the diffusion equation with the spatial boundary conditions is given by 


oo oo 
u(x,t) = }_ bnun(x,t) = Y basin (nrx/ L)e "De? 
n=1 n=1 

The final solution step is to determine the unknown coefficients b, by satisfying conditions on the 
initial dye concentration. We assume that u(x,0) = f(x), where f(x) is some specific function defined 
on0< x < L. Att = 0, we have 

oo 

f(x) = Y by sin (nrex/L). 

n=1 
We immediately recognize this equation as a Fourier sine series for an odd function f(x) with period 
2L. A Fourier sine series results because of the boundary condition f(0) = 0. From our previous 


solution for the coefficients of a Fourier sine series, we determine that 


2 rh , NTX 
bn = iff f(x) sin dx. 
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152 LECTURE 51. SOLUTION OF THE DIFFUSION EQUATION (FOURIER SERIES) 


Problems for Lecture 51 


1. Determine the solution to the diffusion equation ur = Duxy for the concentration u = u(x,t) with 
closed pipe ends corresponding to the boundary conditions ux(0) = ux(L) = 0, and with general 
initial conditions given by u(x,0) = f(x). Use the results of the previously solved problems. 


Solutions to the Problems 


Lecture 52 


Diffusion equation: example 


View this lecture on YouTube 


Example: Determine the concentration of a dye in a pipe of length L, where the dye is initially concentrated 


uniformly across the center of the pipe with total mass Mo, and the ends of the pipe open onto a large reservoir. 


We solve the diffusion equation with homogeneous Dirichlet boundary conditions, and model the 
initial concentration of the dye by a delta-function centered at x = L/2, that is, u(x,0) = Mod(x — L/2). 


The Fourier sine series coefficients are therefore given by 


2 fL L, . nax 
bn = El Moö(x — 5) sin A 


= = sin (nrt/2) 


2Mp/L ifn =1,5,9,...; 
= 4 -2Mp/L ifn =3,7,11,...; 
0 ifn =2,4,6,.... 


With b, determined, the solution for u(x,t) is given by 
u(x,t) = I LEY" sin (Fr gers 
n=0 
When t > L2/D, the leading-order term in the series is a good approximation and is given by 


—T2Dt/L? 


u(x,t) = an sin (71x /L)e 


The mass of the dye in the pipe is decreasing in time, diffusing into the reservoirs located at both ends. 


The total mass in the pipe at time t can be found from 


M(t) = i u(x,t)dx, 


and when t > L?/D, we have 
M(t) = ln, 
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154 LECTURE 52. DIFFUSION EQUATION: EXAMPLE 


Problems for Lecture 52 


1. Determine the concentration of a dye in a pipe of length L at position x and time t, where the dye 
is initially concentrated uniformly across the center of the pipe with total mass Mo, and the ends of 


the pipe are closed. 


Solutions to the Problems 


Practice quiz: The diffusion equation 

1. Determine the eigenfunctions of T' + ADT = 0 with eigenvalues A, = (ny 

2. If u(x,t) = a9/2+ y An cos (n7tx/L) exp (—n?7?Dt/L?) and u(x,0) = f(x), determine the general 
n=1 


formula for dy. 


3. Suppose that a pipe with open ends has an initial dye concentration centered at one-quarter of the 


pipe’s length. Determine the long-time t >> 12 /D solution for the concentration. 


Solutions to the Practice quiz 


155 


156 LECTURE 52. DIFFUSION EQUATION: EXAMPLE 


Appendix 
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Appendix A 


Complex numbers 


View this lecture on YouTube 
We define the imaginary unit i to be one of the two numbers (the other being —i) that satisfies z? +1 = 


0. Formally, we write i= y—1. A complex number z and its complex conjugate Z are written as 
z=x+iy, z=x-iy, 
where x and y are real numbers. We call x the real part of z and y the imaginary part, and write 
x=Rez, y=Imz. 
A linear combination of z and Z can be used to construct the real and imaginary parts of a number: 


1 2 1 > 
Rez = 5 (z +2), Imz = (22). 


We can add, subtract and multiply complex numbers (using i = —1) to get new complex numbers. 
Division of two complex numbers is simplified by writing z/w = zw/ww. The exponential function 
of a complex number can be determined from a Taylor series. We have 


(i0)? (40), Ge)? | (ie) 


a oe ag 
e =1 | (i0) f 21 T 31 t 21 | 51 
02 4 03 g5 
= ( z T a =) -i (6 ats = = cos 0 + isin 0. 
y i Z=x+ly 
wT | 
= | 
£ | 
r | 
| 
0 | 
| 
x 
Re(z) 


The complex plane. 
The complex number x +iy = z can be represented in the complex plane with Rez as the x-axis 


and Imz as the y-axis (see the figure). Then with x = rcos@ and y = rsin0, we have z = x +iy = 
r(cos@ +isin®) = re”. 
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Appendix B 


Nondimensionalization 


View this lecture on YouTube 

An equation can be made dimensionless through a careful exercise in dimensional analysis. First, the 
base dimensions associated with the equation are determined. These may be one or several of the 
fundamental dimensions such as time, length and mass. When working with a dimensional equation, 
specific units must be chosen; for example, length may be measured in meters, feet, miles, or a more 
interesting unit such as light-years (the distance light travels in one year). But when working with a 
dimensionless equation, choosing specific units is unnecessary. 

There is also a very nice theorem called the Buckingham zt theorem, which states that the number 
of dimensional parameters in an equation can be reduced by the number of fundamental dimensions 
associated with that equation through the process of nondimensionalization. In practice, this can 
map seemingly different problems to the same dimensionless equation, and can also vastly reduce the 
number of parameters that must be explored in a numerical solution. 

To illustrate how to render an equation dimensionless, we can consider the pendulum equation 
given by 

ml6 + clð + mg sin09 = Fo cos wt. 


In this equation, the base dimensions are time (t), length (1) and mass (m). The angle 0 expressed in 
radians is dimensionless, since it is defined as the ratio of two lengths (arc length divided by radius). 
The dimensional parameters are m, I, c, g, Fo, and w, for a total of six. The Buckingham 7 theorem says 
that after nondimensionalizing the equation, there should be no more than 6 — 3 = 3 dimensionless 
parameters. 

Each term in this equation must have the same dimension, and from the first term we see that each 
term must have dimension ml/t2, which is the units of a force. The frictional parameter denoted by c 
must have units m/t. The gravitational acceleration g has units of !/ t. 

If we divide the pendulum equation by ml?, and define the natural frequency of the oscillator 
without damping to be wọ = \/ g/l, then the pendulum equation becomes 


i ; F 
0+ 0 +a sind = 2 cos wt. 
m ml 


The units of the various grouping of parameters are given by 


Ca _ ee o 
ll =1/ [wo] = 1/t, 7 =1#,: lea ist 
Nondimensionalization requires a choice of time scale. There is some freedom here, and the choice 
is usually made by deciding which pair of terms should have a coefficient of unity in the final dimen- 


sionless equation. A standard choice is to equate the coefficient of the inertial term 6 (here already 
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unity) to the coefficient of the restoring force sin 0 (here w). In this case, the ws coefficient of the sin 6 
term can be removed in the final dimensionless equation by defining the dimensionless time T to be 


T = wot, 


and using d/dt = (dt/dt)d/dt = wod/dt in the change of variables, resulting in 


d?9 do. 
T2 + 0 +sin@ = ycos PT, 


where the three dimensionless parameters in this dimensionless equation are given by 


c w Fo 
a= , p= Y= $ 
muy wo mlws 


Another choice for nondimensionalization could be found by equating the coefficients of the frictional 


term Ó and the restoring force. Here, one can first rewrite the pendulum equation as 


Tore 6 = Fo cos wt, 
c lc lc 


and observe that the coefficient of the restoring force can be removed in the final dimensionless equa- 


tion by defining the dimensionless time T to be 


aif 
ee ( Ic ) i 
The dimensionless equation then becomes 
0 e dé Siig = 
qa + qq + sind = yoos BT, 
where here 5 
m lcw F 
-eer 
c mg mg 


This latter dimensionless equation is useful when studying very large damping effects, corresponding 


to the limit « — 0. The resulting governing equation then becomes first-order. 


Appendix C 


Matrices and determinants 


View this lecture on YouTube 
A two-by-two matrix A, with two rows and two columns, can be written as 


The first row has elements 41; and a12; the second row has elements a2; and a». The first column 
has elements 211 and a71; the second column has elements 412 and a». Matrices can be multiplied by 
scalars and added. This is done element-by-element and the straightforward definitions are 


[0 #2) _ kay, kan Aii iia) a bu b2) _ (auth am+bn 
an ar kag, kan azn an bza bz az +b an + bog 


Matrices can also be multiplied. Matrix multiplication does not commute, and two matrices can be 


multiplied only if the number of columns of the matrix on the left equals the number of rows of the 
matrix on the right. One multiplies matrices by going across the rows of the first matrix and down the 
columns of the second matrix. The two-by-two example is given by 


aa al (bu bi2\ _ (anbu + 412b21 aybi + a12b22 
azn a2) \br1 bz Azıbıı + 422b91 Azıbı2 + a22b22 


A system of linear algebraic equations can be written in matrix form. For instance, with aj and b; 


given numbers, and x; unknowns, the (two-by-two) system of equations given by 
011%1 + 442%2 = by, — Ag X1 + a22X2 = bp, 
can be written in matrix form as Ax = b, or explicitly as 
a ay f/x) (br 
a a2) Xx b2 
A unique solution to Ax = b exists only when det A ¥ 0, and it is possible to write this unique solution 
as x = AT lb, where A”! is the inverse matrix satisfying ATA = AA~! = I. Here, I is the identity 


matrix satisfying AI = IA = A. For the two-by-two matrix, the determinant is given by the product of 


the diagonal elements minus the product of the off-diagonal elements: 


det A = 411422 — 412421. 
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Appendix D 


Eigenvalues and eigenvectors 


View this lecture on YouTube 
The eigenvalue problem for an n-by-n matrix A is given by 


Ax = Ax, 


where the scalar A is called the eigenvalue and the n-by-1 column vector x is called the eigenvector. 


Using the identity matrix I, we can rewrite the eigenvalue equation as 
(A — AI)x = 0. 


When A is a two-by-two matrix, then 


el |) ge Aa A a 
0 1 421 an — À 


A solution other than x = 0 of the eigenvalue equation exists provided 
det (A — M) = 0. 

This equation is called the characteristic equation of A, and is given by 

A? — (ay + a2)A + (aan — a12a21) = 0, 


which is usually written as 
A? —TrAA+detA = 0, 


where TrA is the trace of the matrix A, equal to the sum of the diagonal elements. The eigenvalues 
can be real and distinct, complex conjugates, or repeated. 
After determining an eigenvalue, say A = M, the corresponding eigenvector vı can be found by 
solving 
(A = Aıl)vı = 0, 


or explicitly in the two-by-two case, 


411 A 412 on Y _ 0 
ary an AM) \ 024 0 


The second equation represented above is always a mutiple of the first equation, and the eigenvector 


is unique only up to multiplication by a constant. 
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Appendix E 


Partial derivatives 


View this lecture on YouTube 
For a function f = f(x,y) of two variables, the partial derivative of f with respect to x is defined as 


Of _ im Se thy) Sy) 
ox h-0 h í 


and similarly for the partial derivative of f with respect to y. To take the partial derivative of f with 
respect to x, say, take the derivative of f with respect to x holding y fixed. As an example, consider 


f(xy) =28y +y°. 


We have of of 
OJ eyz OF _ 443 2 
Sy 6x*y”, T 4x°y + 3y^. 


Second derivatives are defined as the derivatives of the first derivatives, so we have 


and the mixed second partial derivatives are independent of the order in which the derivatives are 


taken: 


ef 2, Df 2 
TE 2x4, 12x°y. 


Partial derivatives are necessary for applying the chain rule. Consider df = f(x+dx,y +dy) — f(x,y). 
We can write df as 


df = [f(x Han, y + dy) — fu y + dy)] + [Fle + dy) FR] = Ldr Lay. 
If one has f = f(x(t),y(t)), say, then 


df _afdx | af dy 
dt dxdt ` ody dt 


And if one has f = f(x(r,0), y(r,0)), say, then 


af ofax ofdy af dafax , af day 
dr oxar' dyar’ 06  ƏxƏð dyad 
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Appendix F 


Table of Laplace transforms 
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APPENDIX F. TABLE OF LAPLACE TRANSFORMS 


5. te” 

6a. sin bt 
6b. sinh bt 
7a. cos bt 
7b. cosh bt 
8. e sin bt 
9. e" cos bt 


10. tsin bt 


11. tcosbt 


12. uc(t) 
13. uc(t)f(t— c) 
14. 6(t —c) 

15. x(t) 


16. ¥(t) 


Table of Laplace Transforms 


Appendix G 


Problem and practice quiz solutions 
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172 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 


Solutions to the Problems for Lecture 1 


a) third order, ordinary, nonlinear 
b) second order, ordinary, nonlinear 
c) second order, partial, nonlinear 
d) second order, ordinary, linear 


e) second order, partial, linear 


Solutions to the Problems for Lecture 2 
1. The modified Euler method is 
1 
kı = Axf(Xn,Yn), k = Axf (Xn + Ax, Yn +k), Yn+1 = Yn + 5 (ka + k2). 


The midpoint method is 


1 1 
ky = Axf (Xn, yn), ko = Axf (xn + ¿AX Yn + 3H), Yn+1 = Yn + k2. 
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Solutions to the Problems for Lecture 3 


2_ 
ij a 


a 
y Y x+4 


b) e” cos y dy = e* (1 + x) dx 


y+1 x 
dy = d 
y Terj” 


c) 


Solutions to the Problems for Lecture 4 


1. Separate the variables, integrate, and solve for the dependent variable. 


Y dy 12 a). 12 2 212 
Pf xdx; 2y =; y Ll goad. 
Ly 1 0 
x dx 1 1 1 = 
b — — dt. =-4 . [Inx- = = 
) Pe = Use in de Inx-In(1-x)] 7 
x(1— x0) _ E xo 


xo(1-x) = xo+(1-xo)et 
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Solutions to the Practice quiz: Separable first-order odes 

1. ay xt2y1/2, y(1) =0. 
dx í 


Y dy PAE 1/2 _ 27,3/2 o er 
han): dx; 2y =3(x =D y= _. 


dy 


Re. _ 
2x Y, y(1) =1. 
y dy x dx. 1 „4 u : 1 
al Bod, er o 1) =1nx; o In x; u re 
dy _ ; _ 
3. = —(sinx)y, y(r/2) =1. 


Yd x : 
f “= -f sinxdx; lny= cosx; y= e®*, 
1 7/2 
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Solutions to the Problems for Lecture 5 


1. 
dy sinx 
a dx x x” 
dy 
) gy JTTX 
2. Letz = 1/x. Then x = 1/z and dx/dt = —(1/22)dz/dt. The differential equation becomes 


—(1/z?)dz/dt = (1/z)(1 — (1/z)), and after multiplying by —z?, dz/dt = —z +1. In standard form, 
the linear ode is dz/dt +z = 1. 


178 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 


Solutions to the Problems for Lecture 6 


a) —=x-y, y(0)=-1. 
dy/dxt+y=x; p= eh% = e; 
y=e*(-1+ fy xe*dx) = e™*(—1 + xe* | — fo etdx) = e7™*(—1 + xe* — e* +1) = x —1. 


2 


dy/dx +2xy =2x; u= elo 2xdx = px”, 


2 
y=e* Le 2xe dx =e Jo «du = eve -1)=1-e*. 


Solutions to the Practice quiz: Linear first-order odes 


1. (14xy +2xy =2x, y(0)=0, 
y2 


1+x2 


x 
K ' x)y] = 2x, (1+ x)|) =f 2xdx = x?, (1 +x? y = x°, y 


Le +2xy =1, y(1)=2, 


d 2 2. |x 2 x a 1+x 
q Py] = Pull = y- f dx =x—1,x°y=1l+x,y= nos 


3.y4+2y=a, y(0)=0, 


u=, y=e™ de ae "dx = — (e* 1) = > (1 ae 


180 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 


Solutions to the Problems for Lecture 7 


k 
1. Use S(t) = Spe + z(a —e™™), with t = T = 40 years, So = 0, r = 0.06, and S(40) = 1,000,000. 


rS(40) _ 0.06 x 1,000,000 
eT-1 e0.06x40 _ 1 
(40 years) x ($5,986 per year) = $240,000, and the amount passively earned from the investment is 


Solving for k, we have k = = $5,986 peryear. The total amount saved is 


$760,000, about three times as much that was saved. 


181 
Solutions to the Problems for Lecture 8 


1. Use v(t) = -8a — ekt/m), where mg/k = 200 km/hr and m = 100 kg. With 

g = (9.8 m/s?) (10°? km/m)(60 s/min)?(60 min/hr)? = 127,008 km /hr?, we find k = 63,504 kg/hr. 
One-half of the terminal velocity for free-fall (100km/hr) is therefore attained when (1 —e7*/™) = 
1/2, ort = mIn2/k = 4sec. Approximately 95% of the terminal velocity (190 km/hr ) is attained after 


17 sec. 


182 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 


Solutions to the Problems for Lecture 9 


1. The current i = dg/dt, with q = CVc. During charging, Vc(t) = E€ (1 =p = and during 
discharging Vc(t) = Ee”'/RC. With i = CdVc/dt, we have during charging, i = (E/R)e'/RC, and 
during discharging, i = —(E/R)e~'/®©. The currents strengths are the same during both processes 


but are in opposite directions. 
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Solutions to the Practice quiz: Applications 


k 
1. The relevant formula is S(T) = U — e "T), where S(T) is the amount saved after T years, r is 


T —rT 
the annual return, and k is the annual deposit rate. Solving for k, we have k = Substituting 


r = 0.1, S(T) = 1,000,000, T = 40, we find k =1,865.74. The total saved is 40k = 74,629.44 ~ $75,000. 


2. Use v(t) = -Ea — ekt/m), where mg/k = 200 km/hr, m = 100kg and v(t) = 150km/hr. With 
g = 127,008 km/hr?, we find k = 63,504 kg/hr. Three-quarter of the terminal velocity for free-fall (150 
km/hr) is therefore attained when (1 — e7™%/") = 3/4, or t = mln 4/k = 7.86 sec © 8 sec. 


3. Use Ve = € (1 = bn. Solve for t to find t = —RC In (1 — E). With R = 3000 Q, C = 0.001 F, 
Ve/E = 0.95, we find t = 8.99s + 9s. 
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Solutions to the Problems for Lecture 10 


1.0=4u ü=-u/g-sind+ fcoswt 


kı = Atf (tn, Xn, Yn), l = Atg (tn, Xn, Yn) 
kz = Atf (tn + At, xn + ki, Yn +h), Ip = Atg (tn + At, Xn + ki, Yn +11) 


1 1 
Xn+1 = Xn + > (kı +k), Yn+1 = Yn + z(h + In). 
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Solutions to the Problems for Lecture 11 


1. With x = xy (t) + xp(t), we compute 


EH p(t)x +q(t)x = (En + ëp) +p (Xp + Xp) +4 (Xp + Xp) 
= (Ën + PXp + Xp) + (Xp + pXp + qxp) 
=0+g(t) 
= g(t), 


since x, and x, were assumed to be solutions of the homoogeneous and inhomogeneous equations, 


respectively. Therefore, their sum is a solution of the inhomogeneous equation. 


186 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 
Solutions to the Problems for Lecture 12 
1. With X; (t) = exp (at) and X2(t) = exp (Bf), we have 

W = (exp (at))(Bexp (Bt)) — (a exp (at)) (exp (Bt)) = (B — a) exp ((a + B)t), 


which is nonzero provided a # B. 


Solutions to the Problems for Lecture 13 


a) r? —1 =0, so that r = +1; 


b) r? +1 = 0, so that r = +i; 


c) 1? — 2r +1 =0, so that r =1. 
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188 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 


Solutions to the Problems for Lecture 14 


1. The characteristic equation is 1? + 4r +3 = (r+1)(r+3) = 0, with roots rı = —1 and rz = —3. The 


general solution is x(t) = cıe + coe *!. The derivative is x(t) = —cje ' — 3c2e7%. Initial conditions 

are satisfied by solving cı + cz = 1 and —c1 — 3c2 = 0. Solution is cı = 3/2 and ca = —1/2. The final 
1 1 

solution is x(t) = 5° : ze St — z ta ze). 

2. The characteristic equation is r? — 1 = 0, with roots rı = 1 and r2 = —1. The general solution 


is x(t) = cıe' + cae *. The derivative is x(t) = cıe' — ce‘. Initial conditions are satisfied by solving 
cı+c2 = Xy and cy — cz = ug. Solution is cy = (xo + uo)/2 and cz = (xg — uo)/2. The final solution 
can be written as 
x(t) = xo cosh t + up sinht, 

where i ; ; ; 

e+e” i e — e 
, sinht = 
2 2 


The function cosht is the hyperbolic cosine of t, and is pronounced ‘cosh tee”, where the “o” sounds 


cosh t = 


like the o in open. The function sinh t is the hyperbolic sine of t and is pronounced ‘sinsh tee’, where 
the i sounds like the i in inch. 
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Solutions to the Problems for Lecture 16 


1. The characteristic equation is r? — 2r +5 = 0, with roots rı = 1+2i and r2 = 71. The general 


solution is 
x(t) =e' (Acos2t+Bsin2t), x(t) = œ (Acos2t + Bsin2t) + 2e' (—Asin2t + Bcos2t). 


The initial conditions x(0) = 1 yields A = 1 and x(0) = 0 yields A+2B = 0, or B = —1/2. The 
solution is i 
x(t) = é (cos 2t — 5 sin2t) ; 


2. The characteristic equation is rT? +1 =0, with roots rı = 1 and ra = -i. The general solution is 
x(t) = Acost + Bsint. The derivative is x(t) = —Asint+ Bcost. Initial conditions are satisfied by 
A = xq and B = ug. The final solution is x(t) = xy cost + uosint. 
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Solutions to the Problems for Lecture 18 


1. The characteristic equation is r? — 2r +1 = (r — 1)? = 0, with repeated root r = 1. The general 
solution is x(t) = (c1 +c2t)e’. The derivative is x(t) = (cy +c2(1 ++4))e*. Initial conditions are satisfied 


by solving cy = 1 and cı + c2 = 0. Solution is cı = 1 and cp = —1. Final solution is x(t) = (1-t)ef. 
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Solutions to the Practice quiz: Homogeneous equations 


1. The characteristic equation is r? — 3r +2 = 0, with roots rı =2,r2 = 1. The general solution is 


21 + cet, with derivative x(t) = 2c1e™% + cze". Initial conditions are satisfied by solving 


x(t) = cıe 
x(0) =cy +c2 = 1 and x(0) =2cı + c2 = 0 to find cı = —1, c2 = 2. The final solution is 


x(t) = —e” 420 = —e” (1-28. 


2. The characteristic equation is r -2r +2 = 0, with complex roots r+ = 1 + i. The general solution is 
x(t) = el (A cost + Bsint), with derivative x(t) = e' ((A+B)cost+ (B — A)sint). Initial conditions 
are satisfied by A = 1 and A +B = 0, or B = —1. The final solution is x(t) = el (cost — sin t). 


3. The characteristic equation is 1? +2r +1 = 0, with repeated root r = —1. The general solution is 
x(t) =e ' (cy + cat), with derivative x(t) = e~' (cy — c1 — cat). Initial conditions are satisfied by cı = 0 
and cz — c1 = 1, or cz = 1. The final solution is x(t) = te. 


192 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 


Solutions to the Problems for Lecture 19 


1. Substituting x(t) = xy (t) + xp, (t) + Xp, (t) into the inhomogeneous ode, we obtain 


d2 
K+ př +qx = gan + Xp, 


d 
Xp) H Pa (xn F Xp, 4 Xp) H q(Xh F Xp, 4 Xp) 


= (Ën + PXp + qXp) 4 (Xp, F PXp, 4 qXp;) (Xp, F PXp, A Xp») 
=0+ 81 +82 
= 81 +82, 


so that the sum of the homogeneous and the two particular solutions solve the ode, and the two free 


constants in x, can be used to satisfy the two initial conditions. 
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Solutions to the Problems for Lecture 20 


1. The homogeneous equation has characteristic equation r? + 5r + 6 = (r+2)(r +3) = 0, with roots 


3. To find the particular solution, we try the 


ri = —2 and rn = —3. Therefore, x, = cie? + ce” 
ansatz x = Ae™. The resulting equation is A — 5A + 6A = 1, with solution A = 1/2. The general 


solution is therefore 


1 1 
x=cıe ee coe cae z a => Me = 302079! = ze 
Satisfying the inititial conditions results in the two equations cı + c2 = —1/2 and —2cı — 3c2 = 1/2, 
with solution c4 = —1 and cp = 1/2. The final solution is 


1 1 1 
x(t) = xe ee ze” = se (1 2e™ 4 ee 
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Solutions to the Problems for Lecture 22 


a) Try x(t) = Acost+ Bsint. Equating the coefficients of cosine and sine, we obtain 5A +3B = —2 
and 3A —5B =0. The solution is A = —5/17 and B = —3/17. 


b) Solve 2 — 32 — 4z = 2e". Try z(t) = Ce’. Determine C = (—5 + 3i)/17 and find 


1 1 
Xp = Re{zy} = Re { 17(-5+33)(cost +isine) | = 77 cost + 3sin t). 
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Solutions to the Problems for Lecture 23 


1. Try x = At + B. Equating the coefficients of tł! and t°, we obtain A = 1 and A +B = 0, or B = -1. 
The solution is xy = t — 1. 
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Solutions to the Problems for Lecture 24 


1. The homogeneous equation has characteristic equation 1? + 3r + 2 = (r+1)(r +2) = 0, with roots 


rı = —land r = —2. Therefore, xp = ce”! + c2e”**. To find the particular solution, we try the ansatz 
x = Ate”?! The resulting equation yields A = —1. The general solution is therefore 
x=ce'+(m—te, ==! + (1-20 +20 *. 


Satisfying the initial conditions results in the two equations cy + cz = 0 and —c¡ — 2c2 — 1 = 0, with 


solution cj = 1 and cp = —1. The final solution is 


x(t) = e™t =the". 
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Solutions to the Practice quiz: Inhomogeneous equations 


1. Try x» = Ae”. Equating the coefficients of the exponential function, we obtain 4A +6A +2A =2, 
y Xp q 8 pP 


1 
with solution A = 1/6. The particular solution is given by Xp(t) = de 


2. Solve the complex differential equation 2 — 2 — 2z = 2e?*, and take xp = Re{zp}. Try zp = ce". 
=3 
Equating the coefficients of the exponential function, we obtain—4C — 2iC — 2C = 2, or C = 1 = A 


We find xp(t =A he —3 + i) (cos 2t + isin 2t) == ren 
p 10 10 


3. Try Xp = At + B to obtain —3A + 2At +2B =t +1. Equating coefficients of powers of t, we find 
1 1 
2A = 1 and —3A + 2B = 1, or A = > B= . The particular solution is given by Xp(t) = zt + 2 
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Solutions to the Problems for Lecture 25 


1. Write the RLC circuit equation as 
L d?g ‚ag 1 Eo 


Rie a eo gO” 


To remove the coefficient of the restoring force in the final dimensionless equation, we define the 
dimensionless time T to be 
t=t/RC. 


After multiplication by RC, the RLC circuit equation becomes 


L da dq 
rg= RCT. 
Rent EoC cos wRCT 


To remove the coefficient of the inhomogeneous term, we can define the dimensionless charge Q to be 


= 4 
Q EC’ 
and the dimensionless equation becomes 
¿PO „fQ +Q = cos pt 
dt? dr ' 
with 
L 
a b = wRC. 


RC 
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Solutions to the Problems for Lecture 26 


1. Write the mass on a spring equation as 


mdx dx k Fo 
PTE } Fr 7*5 z cos wt. 


To remove the coefficient of the restoring force in the final dimensionless equation, we define the 


dimensionless time T to be 
T = kt/c. 


After multiplication by c/k, the mass on a spring equation becomes 


mk d?x _ ax Fo WC 
c2 dt? dr 


To remove the coefficient of the inhomogeneous term, we can define the dimensionless position X to 
be 


kx 
= E 
and the dimensionless equation becomes 
LR IN os BT 
dt? dt das 
with 
mk we 
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Solutions to the Problems for Lecture 28 

1. The long-time solution corresponds to the particular solution. We consider the complex ode 
Z+az+Z= Pt 

with x, = Re{z,}. With the ansatz zp = Ae'P!, we have -B?A+iaßA+A=1, or 


ac l _ (=?) ~ inp 
(=P) Fiap ~ (1 BP apr 


To determine the amplitude of the resulting oscillation, we make use of the polar form of a complex 


number and write 


(1-6) imp = y0- BR + 02 et, 


where tang = -aß/ (1 — B?). The long-time solution is given by 


PERA, 
Xp = Re A — eee e m 


1 
= Nero; cos (Bt + @); 


and the amplitude of the resulting oscillation is therefore given by 


1 
VO = BP + 06 


amplitude = 
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Solutions to the Problems for Lecture 29 


1. F(s) = f e”°' sin bt dt. Integrate by parts letting u = sin bt and dv = e~“‘dt to obtain 
0 


F(s) = E I e “cos bt dt. Integrate by parts a second time letting u = bcos bt and dv = (1/s)e “dt 


b bP b 
to obtain F(s) = z wi ls). Solving for F(s), we obtain F(s) = ane 
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Solutions to the Problems for Lecture 30 


1. The Laplace transform of the ode is given by 


a (2x SXo uo) Hb (sX — xp) +cX =G. 
Solving for X = X (s), we obtain 


G(s) + (as + b)xo + aug 
as? +bs +c ` 


X(s) = 
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Solutions to the Problems for Lecture 31 


1. The Laplace transform of the ode is 


1 


2 
X+55X+6X = , 
u s+1 


with solution 1 


aD (s+1)(s+2)(s+3) 


A partial fraction expansion of X = X(s) is written as 


where the unknown coefficients can be found by using the cover-up method. Multiply both sides by 
(s+1) and sets = —1 to find a = 1/2. Multiply both sides by s +2 and sets = —2 to find b = 1. 
Multiply both sides by s + 3 and sets = —3 to find c = 1/2. Proceeding to take the inverse Laplace 
transform using the table in Appendix F, we obtain 
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Solutions to the Practice quiz: The Laplace transform method 


1. From line 9 of the Table of Laplace Transforms, we have a = —1 and b = 7. The Laplace transform 


1 
is then read off the table as X(s) = enra 


2. Use lines 16, 15, and 3 of the Table of Laplace Transforms to take the Laplace transform of the ode: 
1 
5 X +sX — 6X = 41, Solve for X = X(s) to find X(s) = 


(s+1)(s—2)(s +3) 


1 a b c 
; = = . Using th z hod, fi =1/2, 
3. X(s) Creer stl S40 ae Using the cover-up method, we find a = 1/ 
b = —1,c = 1/2. Then using line 3 of the Table of Laplace Transforms, we obtain 
Lif -x , l -st 
e pze, 
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Solutions to the Problems for Lecture 32 


1. The step-down function is given by 


m 
` 


t< c; 
t-u) =} ne 


The step-up, step-down function is given by 


0, t<a; 
ug(t)-—up(t)=4 1 a<t<b; 
0, t>b. 


2. The Laplace transform is 


Ludo fc) = | eufl- c)dt 
= f FE- oat 

I etl ey ai 

ee dl a f@)ar 

=e SF(s), 


where we have changed variables to t = t — c. 


a) f(t) =t- m(t) (1). 
b) F(s) = L{t} — L[ux(t) (t-1)}. The Laplace transforms are found in Appendix F. We have 
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Solutions to the Problems for Lecture 33 


1. For a > 0, we have | f(x)d(ax)dx = a f(x /a)ô(x")dx'" = 1 f0) = T f(x)ö(x)dx. There- 
fore ö(ax) = IR). For a < 0, the limits of integration reverse upon the substitution, and one finds 


me) 


ö(ax) = or). The general result is ö(ax) = 


x 
2. The integral given by I ö(x’ — c)dx' is zero if x < c and one if x > c. This is the definition of 
0.(x). 


x 
3. The derivative of both sides of 0.(x) = J ô(x' —c)dx’ with respect to x, using the fundamen- 


d 
tal theorem of calculus, results in —6.(x) = 6(x—c). This is only a symbolic expression because 
technically you can not take the derivative of a discontinuous function. 


207 


Solutions to the Problems for Lecture 34 


1. We compute x(1) when the Heaviside function is either zero or one. In both cases we find 


so that x = x(t) is continuous at t = 1. 


2. The Laplace transform is 
1 
2 fn _ 9-278 
s*X(s) +X(s) == (1 e je 


with solution for X = X(s) given by 
1- e_ 275 


X(s) = er: 
Definin 
i 1 
Sry 


and using the table in Appendix F, the inverse Laplace transform of X(s) is 
x(t) = f(t) - un {ft 27), 


where f(t) is the inverse Laplace transform of F(s). A partial fraction expansion yields 


1 1 S 


F = = , 
(s) s(s?+1) s s2+1 


and the inverse Laplace transform from the table in Appendix F is f(t) = 1 — cost. We find 
x(t) = 1 — cost — uzz (t) (1 — cos (t—271)), 


or more clearly, 


x(t) = 


1— cost, ift <27; 
0, if t > 27. 
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Solutions to the Problems for Lecture 35 
1. Taking the Laplace transform of the ode using the table in Appendix F, we have 

sX4X=1-e 2, 


with solution for X = X(s) given by 


1— e 275 
en 
Defining i 
F(s) = re. 
s-+1 


the table in Appendix F yields f(t) = sint, and the inverse Laplace transform of X(s) is 
x(t) = f(t) — uzn (t) f(t — 270) = sin t — un, (t) sin(t — 271), 


or more clearly, 


sint, ift < 27; 
x(t) = ; 
0, if t > 27. 
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Solutions to the Practice quiz: Discontinuous and impulsive inhomogeneous terms 


1. Use the negative of a step-up, step-down function between t = 1 and 2, and a step-up, step-down 
function between t = 3 and 4. We have x(t) = —(u1(t) — uz(t)) + (uz(t) — ua(t)) = —u1 (t) + u2(t) y 
u3(t) —ua(t). 


1 —278s 
2. Laplace transform the ode to obtain s?X — s + X = E d = Solve for X = X(s) to obtain 
s (1 — e7?"s) 1 1 s 1_ cae fl s 
X(s) = f . Writ = , so that X(s) = - -e | - - — |. 
Se s?+1 s(s2+1) di s2+1) s s2+1 an 5° s s*+1 


1, if t < 27; 
Inverse Laplace transform to find x(t) = 1 — u27 (t) (1 — cos (t — 27r) ) = 


cost, ift > 27. 


3. Laplace transform the ode to obtain s?X — s + X =1—e 2°. Solve for X = X(s) to obtain 


s 1 e7 27s 
X(s) = 271 . Inverse Laplace transform to find 


241 s2+1 


. : cost+sint, if t < 27; 
x(t) = cost + sin t — uzr (t) sin (t — 27) = 
cost, if t > 27. 
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Solutions to the Problems for Lecture 36 


1. Iry 
y(x) = L anx", 
n=0 
and obtain M R 
Y n(n —1)anx"? — Y” anx" =0. 
n=0 
In the first sum, shift the summation index downward by two, and then combine sums to obtain 
> (o +2)(n + 1)anı2 — ay a" = 0. 
n=0 
We therefore obtain the recursion relation 


An 
n+2)(n+1)’ 


m= 1n=0,L Lis 


Even and odd coefficients decouple, and we obtain two independent sequences starting with first term 
ay or a1. Developing these sequences, we have for the first sequence, 


1 
a, M= 740 a4 = 12.32 = g” 
and so on; and for the second sequence, 
o 1 o 1 1 
y E SI: Oe ge ep le 


and so on. Using the principle of superposition, the general solution is therefore 


x2 ae Pe 
E ta ta = pan (x+ art 5 = 


= ay cosh x + a; sinh x, 


where 
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Solutions to the Problems for Lecture 37 


1. Substituting y = I _o anx" into the differential equation, we have 


Y"+xy — y= L n(n — Dax"? + ye nayx" — > anx" 
n=2 n=1 n=0 
= Y ((n+2)(n + 1)an42 + (n— 1)an)x” =0. 
n=0 
The recursion relation is 
An12 == (n=1) a 


Starting with ao, we find a = a9/2, a4 = —ag/ 24, and ag = ag/240. Starting with a1, we find a3 = 0, 
as = 0, etc. The general solution to order xê is given by 


x2 6 
via) =a (14 a ba A a) Hax. 
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Solutions to the Problems for Lecture 38 
1. The general solution of the Airy’s equation is given by 

y(x) = aoyo(x) + aryi (x). 


Applying the initial conditions and the known values yo(0) = 1, y¿(0) = 0, yı(0) = 0, and y; (0) = 1, 
we have 
y(0)=1=a, y(0)=1=aj. 


Therefore, the solution is 
y(x) = yo(x) + y1(x). 
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Solutions to the Practice quiz: Series solutions 


1. Using the definitions of the hyperbolic functions, we have 
x —x)\2 x ,-x)\2 
cosh? t — sinh? t = eres (e en 


= (e+e +2) (te 2) =5 14 


oo 
2. Substituting y = L anx” into the differential equation, we have 
n=0 


oo oo 


y +x%y = Y n(n-1)apx?4 Y axt? = Y (n+4)(n+3)anrax"t?+ Y ap? = 0. 
n=2 n=0 n=-2 n=0 


00 
Rewrite as 2a, + 6a3x+ }` [(n +4) (n +3)an+4 + 41] x"*2 =0. Setting coefficients of powers of x 
n=0 
equal to zero results in a) = a3 = 0, 0144 = — 


An 


(n+4)(n +3) 
Starting with a4, we find a5 = —a,/20. Since aa = a3 = ag = a7 = 0, the general solution with terms 


xt » 
up to x9 is y(x) = ay (1 12 | = + ay (: 20 | a 


. Starting with ag, we find a4 = —a9/12. 


oo 
3. Substituting y = ) ' a,x" into the differential equation, we have 


n=0 
y" —xy' +y = Yo n(n— Vay"? — Y nax” + Y anx" = Y (6 +2) (n+1Dan+2 + (1— nn)" = 0, 
n=2 Ve n=0 n=0 
n—1 
Th i lation i = . Starti ith ag, we find aa = —ag/2, dg = —ao/ 24. 
e recursion relation is 4,42 RE er arting with ao, we find a2 ag/2, a4 ag/ 
Starting with a1, we find a3 = 0, a5 = 0, etc. The general solution up to terms proportional to x* is 


2 44 
given by y(x) = ao (1 5 57 =) Hax. 
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Solutions to the Problems for Lecture 39 


1. The odes in matrix form are given by 


OE 


and the characteristic equation of the matrix is given by 


The roots of the characteristic equation are 


a+b y(a+b)9?—4ab+4c?  a+b+y(a—b) +4c? 
s 2 = 2 ' 


Since the discrimanant is non-negative, both roots are real (and distinct provided a Æ b and c #0). 
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Solutions to the Problems for Lecture 40 


0 


—1 
1. With A = ( ) , the characteristic equation is A? + A—2 = (A-1)(A+2) = 0, with roots 
; ; 1 1 
A = 1 and Az = —2. The corresponding eigenvectors are vj = i and v? = a) The general 


1 1 
solution is x = cel ( ) + coe * (2) ‚orxı = cıe + me” and xy = —cye! + 2c7e?, 
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Solutions to the Problems for Lecture 41 


1 


—2 
1. With A= E Al the characteristic equation is à? — 2A +3 = 0, with root A = 14+iV2 


2i 
and its complex conjugate. The corresponding eigenvectors are v = ) and its complex conju- 


v2 
gate. The general solution is constructed from the linearly independent solutions Xy = Re {ve} and 
—2sin /2t 2cos V2t = = 
X = Im !ve!!\, andisx =e! | A +B ‚orxı =e! (—2Asin V2t + 2B cos V2t 
‘i { } v2 cos V2t v2 sin v2t l ( ) 
and x2 =e! (v2A cos V2t + V2B sin v2). 
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Solutions to the Practice quiz: Systems of differential equations 


d (x a b\ fx 
1. The system of equations in matrix form is di ) = l ( ) . The characteristic equation of 
x2 c X2 


the matrix is given by A? — (a + d)A + (ad — bc)A = 0. The second-order ode with the same character- 


istic equation is given by * — (a + d)x + (ad — bc)x = 0. 


d I. 2 
2. The system of equations in matrix form is — Ale = Try x = ve‘! to obtain the 
dt X? 2 1 X2 


eigenvalue problem Av = Av. The characteristic equation of A from det (A — AI) = 0 is 4? — 2A — 3 = 
(A-3)(A+1) = 0, with roots A; = 3 and Aa = —1. The eigenvectors are found by solving 


1 
(A — Ajl)v; = 0, and we find that A; = 3 has eigenvector vı = A and Aa = —1 has eigenvector 


1 1 1 
v = () The general solution is x = c1 () et + cy (3) e”! and the components are xı = 


f 3t 


c1e% + met and x2 = ce% — met. 


d —2 1 
3. The system of equations in matrix form is — A a . Try x= ve! to ob- 
dt X2 —1 -2 X2 


tain the eigenvalue problem Av = Av. The characteristic equation of A from det(A — M) = 0 is 
A? +4A +5 = 0, with roots A = —2 +i and its complex conjugate. The eigenvector associated 


1 
with A is found by solving (A — AI)v = 0, and we find that vı = | . |. The general solution is 
i 


x=e 2 (a Re G) a) +B mf () |) =e! (a e ) +B e and the compo- 
1 1 —sint cost 


nents are xı = e7% (A cost + Bsint) and x) = e *!(—Asint + Bcost). 
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Solutions to the Problems for Lecture 42 
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Solutions to the Problems for Lecture 44 


1. The governing equation may be written as 


md? (x\ (2 1 fx 
kdt xo) | 1 1) (x) 
With ansatz x = ve”, the matrix equation becomes Av = Av, with A = mr? /k.The eigenvalues of A are 


found from det(A — AI) = A? +3A +1 =0, or Ay = -(3- V5)/2 and A, = -(3+ V5)/2. The cor- 


T T 
responding eigenvectors are computed to be vı = (1 (15 + 1)/2) and v2 = (1 (v5 - 1)/2) ; 
The approximate values of the angular frequencies and the eigenvectors are 


k 1 
w1 = 052 E, vı = ; 
m 1.62 
k 1 
un = 162 /£, = ‘ 
2 mv” (e 


220 APPENDIX G. PROBLEM AND PRACTICE QUIZ SOLUTIONS 
Solutions to the Practice quiz: Normal modes 


1. Newton’s law and Hooke's law results in m¥, = —2kx1 — k(x, — x2) and mia = —k(x2— xı). In 


d2 [x —3k k x41 
matrix form, this is written as Mn = á 
dt k —k x2 


2 => 4 
2. The governing equations are given by ue (x) = : ) © ) Try x(t) = ve", and obtain 
X2 = 


Av = Av, with A = T. The characteristic equation of A is A? +4A +2 = 0, with solutions A+ £ 


Thak _ | /(2- = es 


v2. The angular frequencies of the normal modes are given b y | 
5 4 5 y m Je 185, /E. 


-3 1 
3. We find the eigenvectors of the matrix A = ( i | with eigenvalues Aj = —2 + V2 and 


1 14/2 


-1+v2 1 012 
or v = (1+ vV2)vıı = 2.4101]. For Aa, we have = 0, or do = (1- 
r 0 = (1+ V2)on 11 2 i En ( 22 = ( 


1 1 
2)v12 & —0.41v12. Therefore, vı = and v? = ; 
v2)012 12 1 & 2 (u) 


-1-v2 1 
Ay = —2 — V2. We solve (A - A;l)vi = 0 for vi. For A, we have v2 ) (x) =0, 
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Solutions to the Problems for Lecture 45 


a) Computing f(x + 2L), we have 


40 


2L 2L 
(a, cos _ ) + bn sin a ) 


Me 


3 
ll 
a 


f(x +2L) 


(an cos (= + 2n7.) + by, sin (= + 2n7) ) 


| 
n|3 
| 
T 
Ma 


n=1 

ag > NTIX . ATX 
= — + An COS + bn sin 

2 a E 1) 
= f(x). 

b) Computing the average value of f(x), we have 
1 L 
FO) = zg | ‚fo 
1 st ao 2 fan ft NTIX bu fh . nrx 
= | BAK pj nmr 

7 za È (a, E PEST Sin r ax 
ao 
2 


Therefore, a9 = 2 (f(x)). 
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Solutions to the Problems for Lecture 46 


1. If f(x) is an odd function, then its Fourier sine series is given by 
= NTIX 
x) = bn sin ——, 
) L n L 


and f(0 z bn sin0 = 0. If f(x) is an even function, then its Fourier cosine series is given by 
n=1 


a = NTTX 
fo = a È An cos =— : 


n=1 


and 
NICX 


foe} 
= =o 14 sin ——. 


Therefore, f'(0) = -Z Y nan sin0 = 0. 
n=1 
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Solutions to the Problems for Lecture 47 
1. The function f(x) is odd and L = 7. The coefficients of the Fourier sine series are given by 
by = Z [" sinnx dx 
== a cos nx|q 


2 
= — (1 —cosnrr) 


nm 

4 1, ifn odd; 
== x 

AT 0, if n even. 


The Fourier sine series is therefore 


47. _ sin3x  sin5x | 
fla) == (sinx Esa) 


3 5 
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Solutions to the Practice quiz: Fourier series 


1. f(0) = > + O = > de L an 
n= n= 


2. The square wave in the figure is an odd function of period 271, with f(x) = 1 for 0 < x < 71. The 


oo 
Fourier sine series is given by f(x) = } bysinnx with 
n=1 


z m 2 de 4, nodd; 
b, = `f f(x)sinnxdx = =| sinnxdx = —“SOS™* | Im 
TÍO TU JO an lo 0, neven. 
4 sin3x sin5x 
Therefore, = — (sinx4 Eau a 
erefore, f(x) = (sin 3 5 ) 


3. This is the function f(x) of Question 2 where f(x) is shifted to the left by 77/2. If we denote this 
shifted function by g(x), then g(x) = f(x + 7/2). We have 


ds = (sin (x 7/2) 4 mort 7 /2)) Ä KEU m With sin (x + 77/2) = cos x, 
sin (3 (x + 7/2) ) = — cos 3x, etc., we have g(x) = (cosx cons or = 
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Solutions to the Problems for Lecture 48 


1. We consider 
du du 


Ot ax? 
and define the dimensionless time as tT = tD / L? and the dimensionless position as s = x/ L. Changing 


variables, the nondimensional diffusion equation becomes 


du u 


at 98’ 


which is an equation without any dimensionless parameters. 
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Solutions to the Problems for Lecture 49 


1. We have 
ux(0,t) = X'(O)T(t)=0, ux(L,t) = X'(L)T(t) = 0. 


Since these boundary conditions are valid for all t, we must have X’(0) = X’(L) = 0. These are called 
homogeneous Neumann boundary conditions. 
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Solutions to the Problems for Lecture 50 
1. If A = 0, the general solution of X” = 0 is given by 
X(x) = A + Bx, 


and X(0) = 0 yields A = 0, and X(L) = 0 yields B = 0. If A < 0, we write A = —p?, and determine 
the general solution of 
x” = wX = 0 


to be 
X(x) = Ae'* + Be P*, 


The boundary condition at x = 0 results in A + B = 0, and the boundary condition at X = L results 
in Aer! + Be=#l = 0. Combining these two equations yields 


A (er — g) =0. 


The only solution for nonzero u is given by A = 0, and since B = —A, there are no nontrivial solutions. 


2. Nontrivial solutions only exist for A > 0. If A = 0, the general solution of X” = 0 is given by 
X(x) = A + Bx, 


and the boundary conditions X’(0) = X’(L) = 0 require B = 0. We can write the eigenfunction 
associated with A = 0 as Xo(x) = 1. If A > 0, we let A = u?, and the general solution of the ode and 


its derivative is given by 
X(x) = Acos ux + Bsin px, X'(x) = -uAsin px + pB cos px. 


The boundary condition X’(0) = 0 requires B = 0 and the boundary condition X’(L) = 0 results in 
sinuL = 0, with solutions given by fn = nrt/L, where n is a nonzero integer. The eigenvalues are 
therefore A, = (nrı/ L)? for n = 1,2,3,... with corresponding eigenfunctions X, (x) = cos (nrıx/L). 
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Solutions to the Practice quiz: Separable partial differential equations 


1. The units on both sides of an equation must agree. We have [w+] = [D][uxx], or ut?! = [DJul”2, 
where u denotes concentration, t denotes time, and / denotes length. Therefore, [D] = 12171. 


2. To solve uy = c2uxx, we substitute u(x,t) = X(x)T(t) to obtain XT” = EXT. Separating, we 
obtain X’/X = T"/c?T = —A, where A is the separation constant. The two differential equations can 
be written as X” + AX = 0 and T” + Ac*T =0. 


3. The differential equation X” + AX = 0 with boundary conditions X(0) = 0 and X'(L) = 0 has 
nontrivial solutions only when A > 0. We let A = u? > 0, and solve X” + y2X = 0 to find X(x) = 
Acosux + Bsinux. Applying the boundary conditions, we have X(0) = 0 = A and X'(L) = 0 = 
(2n — Dr 
7 
(2n x”) ana 


uB cos uL. Since B 4 0, we must have cos uL = 0. The allowed values of u are given by un = 


forn = 1,2,3,.... The eigenvalues and eigenfunctions are therefore given by A, = ( 


(2n — 1)rıx 
2L 


X, =sin ( )‚fern=123.... 
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Solutions to the Problems for Lecture 51 


1. With the eigenvalues A, = (nrr/ LÝ, for n = 0,1,2,3,..., the differential equation for T = T(t) 
becomes 
re (1272D/1?) r= 


which has solution proportional to 
T Pa 
p= ; 


Therefore, with ansatz u(x,t) = X(x)T(t) and eigenvalues A, we conclude that the functions 


Un(x,t) = cos (nnx /L)e "Dr? 

satisfy the diffusion equation and the spatial boundary conditions for every nonnegative integer n. 
The principle of linear superposition for homogeneous linear differential equations then states that 

the general solution to the diffusion equation with the spatial boundary conditions is given by 


u(x,t) = = + 2 AnUn(x,t) = = + 2 An COS (nrex/L)e "PD LA, 
n=1 n=1 
where following the usual convention we have separated out the n = 0 solution. In the final step, we 
assume that u(x,0) = f(x), where f(x) is some specific function defined on 0 < x < L. Att = 0, we 
have 
a oo 
te) a + Y an cos (n7x/L). 


n=1 
We immediately recognize this equation as a Fourier cosine series for an even function f(x) with 
period 2L. A Fourier cosine series results because of the boundary condition f'(0) = 0. From our 


previous solution for the coefficients of a Fourier cosine series, we determine that 


2 ft NTUX 
An = F f(x) cos A. 
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Solutions to the Problems for Lecture 52 


1. We solve the diffusion equation with homogeneous Neumann boundary conditions, and model the 
initial concentration of the dye by a delta-function centered at x = L/2, that is, u(x,0) = Mod(x — L/2). 


The Fourier cosine series coefficients are therefore given by 


n= N Moö(x ) cos dx 


2 
= a cos (n7t/2) 


2Mo/L ifn=0,4,8,...; 
=! _2M)/L ifn=2,6,10,...; 
0 ifn =1,3,5,.... 


With a, determined, the solution for u(x,t) is given by 
Mo , 2Mo Y 2nTTX An? n° Dt 
u(x,t) = L 1 L > 1)" cos ( L ) exw (A . 


As one would expect for a pipe with closed ends, after a sufficiently long time, the mass of the dye 


becomes uniformly distributed across the length of the pipe. 
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Solutions to the Practice quiz: The diffusion equation 
1. The solution of T' + ADT for T = T(t) up to a multiplicative constant is T = exp (—ADt). Replacing 


n?n2Dt 
12 


A by the eigenvalue A, = (nrt/L)?, we obtain the eigenfunction Ty = exp (- 


2. We have u(x,t) =a0/2+ Y ' an cos (nrıx/L) exp (—n?7?Dt/L?) and u(x,0) = f(x). Therefore, 


n=1 


f(x) = a0/2+ Y ancos(n7x/L). This is a Fourier cosine series for f(x) and the coefficients are 
n=1 


2. Te NTCX 
given by an = a f(x) cos (= )ax 


3. The solution for the concentration in a pipe with open ends is given by 
2,72 
n“n“Dt 


= NTTX 2 rh nx 
,t) = Y) basin | — . With u(x,0) = , we have by = — in | —— )dx. 
u(x,t) L sin ( L ew ( 12 ) ith u(x,0) = f(x), we have a f(x) sin ( T ) x 


L 2 oy% EN . ynrıx 2Mg . ¿nr 
Here, f(x) = Mob (x-7) We have bn == f Mos (x i) sin T )dx = sin (ZZ). The 


_ v2Mo 


2 
leading-order term corresponds to n = 1 where we have bı = a sin (7/4) = . Therefore 


u(x,t) = V2Mo sin (7) exp ( er) i 


L I 


